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FOREWORD 


b 

This digital simulation study, Contract NAS2-7395, was performed under 
the auspices of the NASA Ames Research Center. It has determined the 
performance characteristics of an outer planet's atmospheric probe communica- 
tions link including the transmitter modulator, the RF channel, and the 
receiver/bit synchronizer. The study emphasized the acquisition properties 
and bit error performance of the system under conditions of signal attenuation, 
Doppler, Doppler rate, and turbulence effects which were modeled during a 
previous study for a Saturn/Uranus atmospheric probe system definition, 

Contract NAS2-7328. The simulated flight equipment design was formulated by 
the TRW Systems Group under subcontract to McDonnell Douglas Astronautics - 
East. 

The author expresses his appreciation for the assistance provided by 
Dr. Marvin Stone and Mr. Marshall Huang of TRW Systems Group, Redondo Beach, 
who provided input data and consultation services during the course of this 
study. 
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SIMULATION STUDY OF COMMUNICATION LINK 
FOR PIONEER SATURN/URANUS ATMOSPHERIC ENTRY PROBE 

SUMMARY 

A digital simulation is presented for a candidate modulator-demodulator 
(modem) in a modeled atmospheric scintillation environment with Doppler, 
Doppler rate, and signal attenuation typical of the radio link conditions for 
an outer planets atmospheric entry probe. The results indicate that the 
signal acquisition characteristics and the channel error rate are acceptable 
for the system requirements of the radio link. The simulation also outputs 
data for calculating other error statistics and a quantized symbol stream 
from which error correction decoding can be analyzed. The results from these 
data have not yet been evaluated. 

The receiver and bit synchronizer are modeled in the simulation at the 
level of hardware component parameters rather than at the level of loop 
equation parameters, and the report defines the resultant hardware parameters. 
Likewise, the scintillation is modeled as a low-pass Gaussian process with the 
parameters defined in the report. 

The receiver performance is plotted for nominal hardware parameters and 
a range of signal-to-noise ratios both with and without scintillation. The 
performance is also shown for critical receiver parameter variations. 
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INTRODUCTION 


The objective of this simulation is to determine the performance and to 
refine design parameters for a receiver which serves as one terminal in an 
outer planet relay link. The transmitting terminal is a probe entering the 
planetary atmosphere, and the receiving terminal is in a spacecraft in 
hyperbolic trajectory about the planet. As with most planetary space systems, 
the link is power starved, inferring minimal data rates and the use of error 
correction coding. The link is unique in two respects: First, because of the 

low data rate, the initial frequency uncertainty due to Doppler and oscillator 
uncertainty is about 400 times the data bandwidth. Second, atmosphere 
turbulence may cause significant scintillation on the received signal. 

The initial system design (Ref. 1), Saturn/Uranus Atmospheric Entry 
Probe, conceived a noncoherent link to combat the scintillation at a data 
rate of 44 bits per second, which was rate 1/2 coded to 88 symbols per second. 
Preliminary specifications for the transmit and receive terminals were 
prepared in conjunction with this study (Refs. 2 and 3). In response to 
these requirements, a preliminary hardware design was defined in Reference 4. 
Due to the nonanalytic nature of the effect of scintillation on the receiver 
performance, a simulation was deemed necessary to verify the suitability of 
the design. A simulation program was written to allow a quantized data 
stream to be generated on magnetic tape for subsequent decoding and analysis, 
as well as to investigate receiver acquisition and channel error rate 
performance, but only the latter two programs are described here. 
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CONFIGURATION 


The configuration details described below are the results of the TRW 
Systems Group design from Reference 4, and as amended in subsequent meetings. 
The configuration consists of three pieces of equipment: the modulator, 

the receiver, and the bit synchronizer. 

The modulator, shown in Figure 1, is a continuous phase frequency-shift- 
keyed modulator. This represents a change from the modulator of the initial 
study which was a discontinuous phase frequency-shift-keyed modulator, as pre- 
liminary studies indicated a severe degradation in performance of a discontin- 
uous phase modulator when operating in conjunction with the candidate receiver. 
The data symbols drive a 20 MHz VCXO which is upconverted by mixing with a 
fixed local oscillator to the output frequency. The oscillator will typically 
switch in 20 microseconds; hence, at an 88 symbol per second rate, the modu- 
lation may be considered as true continuous phase frequency shift keying. 

Two other continuous phase modulators were also considered, but discarded 
because of hardware complexity: a three step up-converter and an indirect FM 

modulator. The modulator candidates are discussed in Appendix X. 

The receiver is shown in Figure 2. The receiver signal is first 
converted to an intermediate frequency and the narrow band filtered and 
output to two loops: an automatic gain control (AGC) loop, and an automatic 

frequency control (AFC) loop. The input to the AGC loop drives the coherent 
amplitude detector via an amplifier which is a soft limiter, i.e., it only 
limits when the signal is very large (typically 13 dB above nominal). This 
is to protect the coherent amplitude detector during acquisition, when the gain 
is wide open, at high signal levels. The AGC loop is completed by low-pass 
filtering and signal conditioning to the variable gain IF amplifiers. 

In the AFC loop, the amplifier before the phase detector originally was 
a hard limiter, i.e., its output was always +A. However, early in this study 
it was felt that the AGC should adequately maintain the signal level and a 
soft limiter should be used in its place, performing a function similar to 
that of the coherent amplitude detector amplifier. The AFC loop which tracks 
the incoming frequency is completed by the loop filter and voltage controlled 
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oscillator. 

If the A6C signal is below threshold, a sweep is activated to acquire 
the signal. The input to the sweep from the AFC loop is a polarity sense; 
that is, it reverses the sign of the sweep voltage if the sensed sweep voltage 
exceeds the design level. The design used in the simulation is different 
from what is now proposed in the actual hardware. As used, if the voltage 
to the detector drops below threshold, sweep is immediately initiated. 

Because of the nature of the scintillation, it was decided during the study to 
incorporate a resweep initiation delay into the hardware. Thus, if the 
scintillation caused only a short drop out, the receiver would not start a 
long (in terms of symbols) search operation. The minimum value for the delay 
is determined by the fading time constants. Since the amplitude scintillation 
has a low-pass characteristic with a 3 dB cutoff at 2 Hz (or 4n radians/sec- 
ond), the delay should be greater than 80 msec. The maximum value of the 
delay is constrained by the Doppler rate and the actual loop parameters, 
e.g., the amount of drift the VCO will undergo during the signal absence 
condition. The Doppler rate is 13.3 Hz/sec. Assuming a nominal VCO scale 
factor of 4 KHz/V, phase detector offsets of 1.2 mv, and a loop filter time 
constant t-j of 0.45, the phase detector offset voltage will charge up the 
loop filter capacitor at the rate of aV/t-| =2.7 mv/sec. This means the VCO 
will then drift at approximately 10.7 Hz/sec so the total loop drift, 
including Doppler rates, is of the order of 24 Hz/sec. Since the receiver has 
a lock-in range of 38 Hz (af L = 2^f n ), the delay should be less than 1.6 sec. 
Given these two extrema, a time constant of 0.5 sec is recommended to delay 
the start of sweep after the lock indicator indicates loss of lock. This 
value could be varied if tests indicate a better value. 

Mixer offsets were not simulated since the mixer output drift was only 
+0.15% over the anticipated voltage range. 

A square foot function was considered proper and implementable for the AGC 
signal conditioner, since the proposed device was linear over a 50 dB range, 
and the receiver requirements only covered 20 dB, i.e., nonlinearities in the 
implementation device need not be simulated. 
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Figure 3 shows the bit synchronizer. This is a conventional early-late 
gate synchronizer with the addition of a baseline correction circuit to 
"correct 0 the effective drift of large numbers of identical bits. As in 
the case of the receiver proper, drift of the multipliers was determined 
to have only a 4% timing error, and was judged to be small enough to be 
idealized as zero for the simulation. A 1-Hz loop noise bandwidth was chosen 
to maximize the signal-to-noise ratio in the loop while maintaining reasonabl 
acquisition times. 



SOFTWARE 


In order to simulate the hardware as closely as possible, the approach 
taken is to simulate the lowest feasible component. For example, each 
multiplier and filter in a loop are individually simulated, rather than an 
expression describing the operation of the loop. 

A powerful simulation tool is the concept of complex amplitude notation. 
In this concept, used in the simulation, and described more fully in Appendix 
I, the "carrier" is normalized out of the simulation. Generically the signal 
is of the form 

x(t) = A cos (w Q t + e), 

wherein A is the amplitude, w Q the carrier and e the phase. By suitable 
normalization, the signal is thus represented in complex notation. as 

x(t) = A exp (je). 

As an example, a frequency-shift-keyed signal (used herein) has two tones, 

x-j (t) = A exp ( j ttA f t ) 
and 

X 2 (t) = A exp ( - jrrAft) 

where Af is the frequency separation between tones. 

Amplitude and phase noise are simply introduced as 

x(t) = A(1 + N ft ) exp [j (e + N p )]. 

One of the principal reasons for this simulation was the determination 
of the effects, on the modem, of atmospheric scintillation. Atmospheric 
scintillation arises from localized variations in the refractive index of 
the atmosphere, chiefly occurring due to temperature variations. Given that 
a "blob" of atmosphere occurs, it successively breaks into smaller and smaller 
"blobs" until its energy is dissipated. With numerous blobs across the 
propagation path, the electromagnetic ray is diffracted in a random fashion, 
resulting i'n variations in the signal amplitude and phase. This simulation 
has modeled amplitude and phase scintillation as 

A(t) = A q (1 + x(t) ) 
and 

e(t) = 9 q + e-|t + e 2 t 2 + s(t), 
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i.e., a linear variation of the amplitude with its scintillation, x(t), and 
the phase with its scintillation, s(t). The e^t and e^t^ terms represent 
Doppler and Doppler rate, respectively. The x(t) and s(t) terms are taken 
as low-pass zero-mean Gaussian processes with root variances of 0.23 and 
0.47, respectively. The filter corner frequencies are 2. and 0.2 Hz, 
respectively and both roll off in proportion to frequency to the -4/3 power. 

The 4/3 f rolloff presents significant modeling effort, and is detailed 
in Appendix II. It was seen that the "more efficient" Z transform approach 
yielded little or no benefit for a 4/3 order filter, in marked contrast to 
its effect on integer order filters. Thus, tapped delay lines were 
employed in the computer routine. To decrease computer time, the simulation 
linearly interpolates the scintillation between "tap times". 

Figure 4 is a block diagram of the simulation receiver. This is "identi- 
cal" to the hardware block diagram with two exceptions: first, the 

simulation receiver starts with the input to the first mixer, the radio 
frequency amplifier being only a constant; and second, the soft limiters 
before the coherent amplitude detector and phase detector are not modeled as 
they only serve in the hardware to protect the following devices from 
excessive power. It may be argued that the soft limiters would have an effect 
during high signal -to-noise acquisitions. However, this region of operation 
was judged to be inappropriate in this study. It should also be noted that 
the AGC signal processor, a square root device, is treated specially. In 
the routine, because the square root function is one of the more time 
consuming functions, the square root is linearly interpolated between 
logarithmic values yielding a maximal 1.46% error. For similar reasons the 
trignometric functions are rounded to one degree and tabulated. Note also 
in Figure 4 that a soft output is added to aid in coding analyses. A 
rationale for setting the quanta levels is found in Appendix III. 

An abbreviated flow chart for the error rate routine is presented in 
Figure 5. After the initial conditions are established (data, dimension and 
complex statements) the conditions for this problem are read in. Typically 
these are variations in signal strength, bandwidths, etc. Next the loop 
conditions are calculated (time constants, etc) and the turbulence delay lines 
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initialized. The routine then generates samples (typically 40) for each bit 
and the loops (AGC, AFC, and bit synchronizer) exercised. The receiver 
starts below threshold, i.e., the AGC/AFC must acquire and the nimber of bits 
to initial lock is determined. The error counter begins when the receiver 
and bit synchronizer have locked, typically after the reciprocal of the 
natural frequency of the bit synchronizer loop or approximately 100 bits. 

This eliminates the error rate "contamination" by both receiver and bit 
synchronizer lockup effects. Additionally, three other statistics are 
tabulated: the bit spacing between errors; the conditional probabilities of 

each quanta level, given the "preceding" level of the soft decision; and the 
number of times the AGC signal drops below threshold and the number of bits 
in that dropout. The printout gives the average number of bits per dropout. 
The routine then continues until the specified number of bits has been pro^ 
cessed. A complete list of symbols is found in Appendix IV, and a complete 
listing in Appendix V. The language is standard Fortran IV. No special 
computer library subroutines are required. 

A typical error rate printout is found in Figure 6. Block 1 has the 
inputs (signal level, bandwidths, etc,). Block 2 contains intermediate 
information such as VCO gains and time constants. Block 3 has the mean error 
rate data: line 1 for the (baseline) filter and sample, and line 2 for an 

integrate and dump output. For this run of 9.4 dB E/No, the mean error was 
0.0247 with a standard deviation of 0.0035. The mean absolute sample level 
was 0.456 with a standard deviation of 0.205. Reading in block 4, lockup 
occurred in one bit, but there was one AGC dropout four bits long. Block 5 
is the error spacing matrix, where for this run 8 bits were in error 1 bit 
spacing apart, 3 bits 2 spacing, 1 bit 3 spacing ... 1 bit 96 spacing and 
4 bits with 100 or greater spacings. In this run a moderate amount of 
"burstiness" occurred, i.e., errors not randomly spaced. Block 6 is the data 
for four conditional probability matricies, which are read by the figure in 
the block. The columns of the matrix signify the quantized level of the "B" 
symbol (current symbol) and the rows signify the level of the "A" symbol 
(preceding symbol). The most negative levels occur in the first column and 
the first row, or the upper left position. For example, given the input 
was -1, -1, which occurred 453 times, an error occurred 22 times. For this 
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matrix the -1, -1 quadrant contains all the "corrects", the +1, -1 and -1, +1 
quadrants single errors and the +1, +1 quadrant double errors. Reading in the 
-1 , +1 quadrant when B (current bit) was position 5 (in error) the last bit, 

A, was four times in 1, two times in 2, two times in 3 and two times in 4. 
Finally block 7 has the end points of the random number generators plus a 
monitor recommending the best sampling time (phase of the bit synchronizer in 
sample counts). 

As the complete receiver simulation is somewhat lengthy, plus having 
special steering and routines for error measurement, a special routine was 
adapted for acquisition studies with its own steering requirements. An 
abbreviated flow chart is shown in Figure 7, with details in Appendix VI. The 
initial setup, problem reading and conversions are identical to those of the 
error rate routine. In this case, however, Doppler is injected, and once the 
receiver acquires, another “case" is run to determine the probabilities of 
acquisition under various conditions. It should be noted on the simulation 
block diagram that the Doppler is injected as a loop stress, rather than a 
phase time product at the input. Comparing the loop responses in the two cases 
reveals identical baseband performance even though the double frequency terms 
which are filtered out are slightly different. Returning to the flow chart, 
after each acquisition trial, intermediate results are printed. For each 
trial the receiver is reinitialized and samples are generated. If the sweep 
goes past the input Doppler, a miss is registered. If the AGC signal goes 
above threshold an acquisition may have occurred. A check of whether the 
sweep got to the Doppler determines whether the condition was a true or a 
false acquisition. 

A typical acquisition printout is shown in Figure 8. Blocks 1 and 2 are 
as stated before. Block 3 shows the conditions at the end of each trial. In 
the first 12 sweeps the receiver acquired. In the 13th sweep it missed 
acquisition. From the 14th to the 20th trials it acquired, but missed again on 
the 21st trial. In no case was a false acquisition found. The fifth column 
contains the bit count at the end of that run, the following three columns 
give the sweep voltage, its value after the AFC loop filter, and the value 
input to the bit synchronizer. 
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During the "fabrication" of the acquisition software, it was found that 
serious aliasing could occur if the sample rate were insufficient. Referring 
to the input mixer in Figure 4, it is seen that as far as the computer is 
concerned it is merely multiplying two strings of complex numbers. If they 
appear to correlate after filtering, large correlated numbers (rather than 
random zero mean numbers) are detected in the CAD circuit signifying lock. 

Now, if the samples are widely spaced in the time domain, apparent correlation 
due to aliasing can occur. This happens periodically when uniformly sweeping 
the VCO, resulting in a false acquisition. Because of the relatively low 
sweep rates, this condition occurred long enough for the CAD filter to 
respond. During the error rate studies a sample rate of 40 samples per bit 
was used. When this rate was used in the acquisition studies, false 
acquisitions always occurred. The sample rate had to be increased to 200 
samples per bit for 1 kHz Doppler before the aliasing effect was negligible. 
Higher Doppler shifts would require correspondingly higher sampling rates. 
Therefore, the acquisition studies were run at 1 kHz Hertz Doppler to 
conserve computer time. 

Both of the routines have seven subroutines: 

NRANDM: Gaussian Number Generator 

RANDM: Uniform Nunber Generator 

B1FILT: Bandpass Filter 

BLOBA(P) : Amplitude and Phase Scintillation Lines 

TRIG: Trigonometry ( 1 —deg table for computer speed up) 

ASQRTF: Square Root (linear interpolation for computer speed up) 

DATAS: Data Source (six-bit feedback shift register) 

Some subroutines used in the checkout phases, but not used in the 
production runs, are given in Appendix VII. 
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RESULTS FOR NOMINAL CONDITIONS 


The first step in the production runs was to determine the performance 
for nominal conditions. These are given in Table I. All the performance data 
points are based on 1900 bit scintillation runs. Before this simulation 
could be done, however, two additional parameters had to be found: the proper 

threshold level and the AGC gain. 


TABLE I 

NOMINAL CONDITIONS 


Parameter 

Symbol 

Value 

Bit Rate 

BRATE 

88. bps 

Minimum Loop Signal 

ENOMIN 

9. dB 

Tone Separation 

DELF 

62. Hz 

IF Bandwidth 

BIF 

1000. Hz 

AFC Loop Bandwidth 

BL2 

176. Hz 

AFC Phase Error 

PHIDEG 

10. deg 

AFC Pull-In Range 

FDIFF 

342. Hz 

Bit Sync Bandwidth 

BL2B 

1. Hz 

Bit Sync Phase Error 

PHIB 

10. deg 

Bit Sync Pul 1 -In 

FDLFFB 

5. Hz 

AGC Filter Constant 

TAU3 

.07957 

Sampling Filter Constant 

TAU4 

.001808 

Baseline Filter Constant 

TAU5 

.01808 

Doppler Filter Constant 

TAU6 

.10603 

Baseline Filter Gain 

A5 

.5 


Normally the threshold is initially set at some factor greater than the 
noise-only output, and then varied to minimize dropouts and minimize false 
acquisitions. For the simulation, the probabilities of a voltage occurring 
with noise only (herein -20 dB E/N , since it was easier to reduce the signal 
than to eliminate it) and with signal was determined. This is shown in 
Figure 9. All of the data referred to in this section appear in Appendix 
VIII. It should be noted that these early runs employed a slightly different 
printout format with more internal diagnostics than used in the production 
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runs. The symbols are all identical however; thus no confusion is anticipated 
in reading the data. Based on this information a threshold of 6.5 was used. 

The AGC signal processor attempts to keep the voltage into the phase 
detector constant, to keep a constant open loop gain. This is done by setting 
the IF gain according to the reciprocal of the square root of the filtered co- 
herent amplitude detector output. For a noise-only input, the gain is 
maximum, and a gain function ABC on Figure 10 is produced. WhiTe determining 
this function during the debugging operation, it appeared that additional gain 
below "minimum signal" would be beneficial, particularly in a scintillation 
environment, e.g., curve ADE or curve AFG. The additional gain herein is 
termed the dynamic range, DYNR. Results of such runs are shown in Figure 11. 

A DYNR of 3 dB appears to offer the best performance with scintillation, and 
was thus chosen for the nominal conditions. 

The first parameter of interest is the mean error rate. This is shown in 
Figure 12. The vertical marks about each experimental point indicate the 
measurement standard deviation. Note that these runs count errors even 
during AGC dropout after initial acquisition. The mean error rate without 
scintillation is well within the "specified" value in the receiver specifica- 
tion, Ref. 3. With scintillation the performance is acceptable to about 
10.5 dB. It should be noted that this performance refers only to the mean 
error rate. In a coded system the distribution of errors is also important. 

The data dumps in Appendix VIII show that the nonscintillation error 
distributions are roughly random (the error history matrix), but the 
scintillation runs have significant burst errors, i.e., groupings of errors. 
Further coding analysis is required to evaluate the performance of the 
system with scintillation. 

Further study of the conditional probability matrices should enhance the 
evaluation of the system's capability to support coding. As noted in the pre- 
vious section, this is moderately influenced by the soft quantization level. 

The bulk of the data was generated prior to the coding analysis of Appendix III 
and used a soft quanta of ASOFT/2 = 0.15. Figure 13 shows the mean and 
standard deviation of the sampled level voltages of the data of Figure 12. 

This was diagnostic information to determine the quantization level. At an 
E/N q of 9.4 dB, the nominal minimum, the standard deviation is 0.2; thus. 
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ASO FT/ 2 = 0.15 corresponds to a quanta level of 0.75 (AS0FT/2/std. dev.) of 
the standard deviation. According to Appendix III, 0.6 would be more 
appropriate, or AS0FT/2 = (0.6 std. dev.) or 0.12. 

Another interesting factor seen in Figure 13 is the slight increase in 
"output" signal -to-noise ratio with the input. At E/N Q of 7 dB we have a 
mean to standard deviation of 0.45/0.23 (5.83 dB) and at E/N Q of 11 dB we 
have a mean to standard deviation of 0.47/0.19 (7.87 dB). There is only a 
2.04 dB apparent increase in output for a 4 dB input increase indicating 
some memory, intersymbol interference, or possibly implied "noncoherent" 
signal detection, i.e., autocorrelation of signal plus noise. 

The second parameter of interest is the probability of acquisition. Runs 
were made at 9.4 and 12.4 dB, with and without scintillation. The data in the 
appendix indicate 100% acquisition. At lower signal to noise, as in Figure 8, 
missed acquisitions occur. There was insufficient data, however, to plot the 
function. 
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RESULTS FOR PARAMETER VARIATIONS 


A number of runs were made to investigate the effects of varying the para- 
meters of the receiver. These may be roughly broken into two cases: system 
variations, and hardware variations. Table II lists the system variations. 

The data for this section are given in Appendix IX. 


The results of modulation index variation is shown in Figure 14. 

As 

expected the optimum is 0.7, with 

and without scintillation. 



TABLE II 


SYSTEM 

PARAMETER VARIATIONS 


PARAMETER 

VALUES 


DELF (Modulation index oriented) 

26.4, 44, 62, 88, 132 

Hz 

BIF (IF Bandwidth) 

500, 750, 1000, 1500, 2000 

Hz 

BL2 (AFC 2B l ) 

88, 100, 164, 176, 188 

Hz 

TAU4 (Sampled Output Filter) 

88, 144, 200 

Hz 


The second system parametric, the intermediate frequency bandwidth, is of 
some implementation importance. If the bandwidth is too narrow the filter 
could produce sufficient phase shift in the loop to cause instability. An 
intermediate frequency bandwidth of ten times the loop filter bandwidth, or 
1 kHz was initially recommended. Because of the high Q of this device, wider 
bandwidths might be easier to fabricate. Figure 15 illustrates the data. 

These represent a second set of runs, because an initial set was required to 
establish the proper level at the sampled output and to define the simulation 
sampling requirements in the scintillation. The runs at and below 1 kHz were 
run with 40 samples per bit, and those above proportionally higher. Without 
scintillation, the higher bandwidths give slightly better performance; 
apparently beginning to level out in the vicinity of 1500-2000 Hz. With 
scintillation, the "knee" appears to be between 1000 and 1500 Hz. 

The results with the receiver loop bandwidth parametric are shown in 
Figure 16. The nominal 176 Hz case is optimal both with and without 
scintillation. 

The final system parametric is the variation of the corner frequency 
of the sampling filter. In conjunction with this, the output of the 
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synchronizer integrate-and-dump was also measured. These data, in Figure 17, 
are also the second set of runs, the first being required to optimize the 
sampling time. It is seen that wider bandwidths degrade the performance, 
although less dramatically when scintillation is present. The integrate- 
and-dump output appears to be uniformly better. Incidentally, the three 
points shown for each integrate-and-dump output represent one datum; since no 
parameter change is made for it; the apparent variations are experiment scatter. 

The expected hardware variations are summarized in Table III. Of these, 
only ENODB, DELF, and BL2 are judged significant. Two "worst cases" were 
defined; ENODB of 9.2, DELF of 66, and BL2 of 164; and ENODB of 9.2, DELF 
of 60, and BL2 of 188. The simulation results are shown in Figure 18, 

Worst case 2 shows little variation from the nominal, within one sigma of 
the measurement. Worst case 1 shows between 0.5 and 0.75 dB degradation. 

TABLE III 

HARDWARE VARIATIONS 


Parameter 

Nominal 

Min 

Max 

Comment 

ENODB 

9.4 dB 

9.2 

9.9 

Based on -142.8 dBm, noise at 
-174 dBm and NF of 2.3 dB. 

DELF 

62 Hz 

60 

66 

Oscillator pulling range accuracy 

BIF 

1000 Hz 

- Negligible - 

Crystal manufacturer data 

BL2 

176 Hz 

164 

188 

RSS estimate of gain changes, and 
scale factors 

SVRHZ 

2.2 kHz/sec 

2.1 

2.3 

5% variance to sweep and VC0 
scale factor 

EVMAX-EVMIN 

35 kHz 

33 

37 

5% sweep drift and VC0 scale 
factor 

TAU3 

.07957 

.07718 

.08196 

3 % RSS time constant 

TAU4 

.001808 

.001754 

.001862 

Same as TAU3 

TAU5 

.01808 

.01754 

.01862 

Same as TAU3 
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CONCLUSION 


Digital computer simulation of a data link has shown that the candidate 
modem appears, from the mean error rate viewpoint, to be acceptable to the 
systems requirements. A complete verification of the usage of the modem for 
a coded telemetry system must await further coding analyses. The performance 
is well within the specified error rate bounds without atmospheric scintilla- 
tion, and the errors appear to be random. Performance is also acceptable 
with scintillation but the errors appear to be bursty. 

The system parametric studies reveal that: 

o a slightly larger intermediate frequency bandwidth of 1500-Hz 
would be more advantageous than the current 1000-Hz nominal. 

o the data output sampling filter of the bit synchronizer could 
be improved, e.g., a straight integrate-and-dump, or perhaps a 
higher order filter. 

o the effect of hardware degradations are within acceptable range. 

o in the range of signal strengths anticipated by the Reference 1 
system studies, the receiver appears to acquire every time. 

Additional runs are required to fully trace the probabili ty-of- 
acquisition curve. 

Because of scintillation modeling uncertainties, it appears worthwhile 
to try other scintillation models and also to explore the effects of the 
resulting error distributions in coding studies. 
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FIGURE 1 

SINGLE CONVERSION DIRECT FM MODULATOR 
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FIGURE 2 

RECEIVER DETAILED BLOCK DIAGRAM 
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FIGURE 4 

SIMULATION BLOCK DIAGRAM 
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FIGURE 7 ABBREVIATED ACQUISITION FLOW CHART 
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FIGURE 8 

TYPICAL ACQUISITION PRINTOUT 
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FIGURE 9 

AGC VOLTAGE DISTRIBUTION 
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FIGURE 11 

MEAN ERROR RATE VERSUS 
AGC DYNAMIC RANGE PARAMETRIC 



FIGURE 12 

MEAN ERROR RATE VERSUS E/N 
FOR NOMINAL CONDITIONS 0 
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FIGURE 13 


FILTER AND SAMPLE OUTPUT 
VOLTAGE VERSUS E/N„ 










FIGURE 14 

MEAN ERROR RATE VERSUS 
MODULATION INDEX PARAMETRIC 




FIGURE 15 

MEAN ERROR RATE VERSUS 
B1F PARAMETRIC 
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FIGURE 16 

MEAN ERROR RATE VERSUS 
BL2 PARAMETRIC 


33 




FIGURE 17 

MEAN ERROR RATE VERSUS 
TAU4 PARAMETRIC 
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FIGURE 18 

MEAN ERROR RATE UNDER TWO 
WORST CASE HARDWARE VARIATIONS 
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APPENDIX I 

COMPLEX AMPLITUDE NOTATION 

Complex amplitude notation is a convenient way to normalize the notation 
to a reference frequency. The signal waveform is thus 

x(t) = Re (x(t) exp (j w Q t)}, 

where the reference frequency is w Q , and the complex amplitude is 

x(t) = A exp ( j e ) = A cose + j A sin e. 

Herein, A is the signal amplitude, and e the signal phase. The conventional 
notation is apparent by combining, 

x(t) = Re { A exp (je) exp (j w Q t) > 

= A cos (w t + e). 

When multiplying in complex amplitudes 

z(t) = x(t) y(t) 

or z(t) = Re (x(t) exp(j.w Q t)} Re {y(t) exp (j w Q t) }. 

Using the identity 

Re {x ( t) > Re { y (t) } = |fte { x(t) y(t) + x{t) y*(t)], 
the star denoting the complex conjugate, we have 

z(t) = ~ Re {$(t) y(t) exp (j (w q +w 1 )t) 

+ x(t) y*{t) exp {j (w 0 - Wl )t)}. 

Normally one beats down, so the complex amplitude of the difference frequency 
term is ^ x(t) y*(t) . 

When filtering x(t) with h{t) to yield y(t), one normally expresses this as 
y(t) = x{t) h(t-x)d . 

— oo 

In complex amplitude notation 

y( t ) “ C Re{x(t) exp (j w q t ) } Re (ft{t-x) exp (j w^t-x))} dr 
= ^•Re{f^x( , r) fr(t- T ) exp( j (w Q -w^ )t )exp (j w^t) dx } 

+ ^-R e{/j?(x) fr*(t-x) exp (j(w 0 +w-|) T ) exp(-j w^Jdx }. 
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Normally, x(t) and h(t) are referenced to the same frequency, so that 

y(t ) = j Re {exp (jwt) /“ x(x) ft(t-x) d t) 

+ j Re {exp(-jwt) x(x) ft*(t-x) exp (2 jwt ) d T >. 

When the bandwidths of the signal and the filter are narrow compared to the 
center frequency, one may neglect the second term, so 

y(t) = \s ; S( t) ff(t-r) dt. 

Finally, when the carrier is zero, to represent baseband signals, 

y(t) = J(t) ft(t-x) dx. 

Figure 1-1 illustrates the usage of complex amplitudes in the simulation. 
The input has the signal amplitude M Q and phase e-j plus noise e . This 
represents the radio frequency, and is beat down by the first mixer to inter- 
mediate frequency by straightforward multiplication. After filtering, the 
extraneous phase is represented by <p , with M* representing the filtered 
amplitude. The real part of the signal from the phase detector is just the 
amplitude times the sine of the phase. ^Similarly the real part of the signal 
from the coherent amplitude detector is just the negative of the amplitude 
times the cosine of the phase. 



FIGURE 1-1 

COMPLEX AMPLITUDE USAGE 
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APPENDIX II 

SCINTILLATION FILTERING 


ABSTRACT: A tapped delay line filter is presented to represent atmospheric 

scintillation effects. Data is given to show the aliasing effects 
of the truncation. 


1. Introduction - It is well known that atmospheric turbulence, herein defined 
as variations in the complex dielectric constant, can cause variation in the 
amplitude and phase of a received electromagnetic wave. From Tatarski in 

Ref (A2-1) , the power spectral densities of the amplitude and the phase is 
proportional to the -8/3 power of the frequency above cutoff. This assumes the 
Kolmogorov form of the turbulent flow. The nature of the fluctuations is dif- 
ferent however for amplitude and phase; amplitude fluctuations being dominated 
by the microscale of the turbulence, and phase fluctuations being dominated by 
the macroscale of the turbulence. 

According to Ames Research Center RFP, the atmospheric effects 
can be adequately modeled by assuming the amplitude effects of the form 

A(t) = A q (1 + x (t)) 

where x(t) is a Gaussian process of zero mean and a root variance of 0.23, 
low pass filtered with a cutoff of 2.0 Hertz and a f“^/3 roll off (note 
herein the signal is a voltage term). Similarly, the phase effects could be 
modeled as 

e(t) = 6 + e.t + 9,t 2 + s (t) 

O 1 l 

where s(t) is a Gaussian process of zero mean and a root variance of 0.47, low 
pass filtered with a cutoff of 0.2 Hertz and a -4/3 roll off. 

This note presents the digital form these filters will take. 

2. Filter Form - In electrical engineering terms the transfer function of 
the filter is 

F(s) = (a/(s+a)) 4/3 , 

where a is the cutoff frequency in radians, and s is joj. 

By simple algebra, the comer is then 
3/4 

to/a = (2 ' - 1) = .8257074727, 

or for amplitude a = 15.21891351 
and for phase a = 1*521891351 . 

S 


n-i 



The Laplace transform of the Impulse response from Ref (A2-2), is 

h(t) = [a 4/3 /r(4/3)]t 1/3 exp(-at), 

and T(4/3) - .892979511 . 

For t of zero or infinity of course 

h(t) | = 0 
t = 0 
t = » 

The function has a maxima at 

dh(t)/dt | = 0 

t = JL 
3a 

The impulse response is shown in Figure II-l. 

3. Z Transform - The quickest computational form of a filter is generally the 
Z transform (Kef (A2-3)). This is just the Laplace transform with the substitution: 

z = exp(sT) 

Utilizing the simplest interpolator, zero order hold, the discrete transform is 
(1 - z" 1 ) Z (j F(s)). 

Since for the case in question 

^((a/(s+a)) 4/3 ) = a 4/3 r _1 (4/3) t" 1/3 e~ at , 
then £($>(.)) - / C a 4/3 r -l (4/3) ,-1/3 .-at d , . f(t) 
and by definition 

00 

ZO^-F(s)) - E fCnt)z‘ n 
S n=0 
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Thus 


oo 


Z(^F(s»- 


z 

n=0 



r' 1 (4/3) / t (nT)- 1/ V ant d(nT)}z-“ 
0 


From Abramowitz and Stegun, Ref (A2-2) , the above incomplete gamma function 
can be written as a Chi-Squared function, or 

OO 

Z(|f( 9 )) = Z {P(2ant||)}z~ n , 
n=0 


Then the complete transform is 

CO 

(l-z _1 ) Z (-F(s))= (l-z'bt £ P(2anT||)z" n } 

n-0 


= P(o]|)z"° + P(2aT|-|)z -1 + P(4aT||)z -2 +. . . 

- P(2aT||)z' 2 - P(4aT||)z” 3 -. . . 

= 0 + P(2aT||)z _1 + (P(4aT||) - P(2aT|-|) }z~ 2 
+ {P(6aT[|) - P(4aTj|)}z" 3 
+ . . . 

The output $ for the input I is then 

♦ - P(2aT||)l((n-l)T) + (P(4aT||) - P(2aT||)}l((n-2)T) 

+ {P(6aT |) - P(4aT||)}l((n-3)T) 

+ . . . 

Since there appears to be no more concise form, at least to the writer, 
in this case it appears that the Z transform is no better than a conventional 
delay line representation, 
as a series, 

p(x 2 N) = r _ 1 (u/2) 

the conventional delay line 


Also, since the Chi-Squared function is evaluated 

j 7+ n 

" (-D (X 2 /2) 2 

n = 0 , ,u , 

n 1 ( j + n 

representation appears preferable. 
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4. Delay Line - There are two design parameters which must be set when representing 
a filter as a tapped delay line: the tap spacing, T, and the line length, NT. 

The former describes the "graininess" of the impulse response, and effectively 
defines the highest frequency represented, 

f = a/(2y) = Kf./y 
m d 

Therein T = y/a, 

f^ is the corner frequency, and for the 4/3 filter 
K = 3. 60472837S. 


The latter describes the fold over aliasing, and effectively defines the lowest 
frequency (excluding zero) represented, 

f L = l/'CNT) 

As a rule of thumb, the maximum tap contribution, Th(t) , should be of the order 
of 0.1, or 

T * <0.1)3 1/3 r(4/3)/(ae _l/3 ) 

» .179740828/a . 


Similarly, the “length" should sura to 0.95 to 0.99, i.e. , 
question 


P(2aNT||) = a 4/3 T 1 (4/3) f 


t sa 2aNT 


<k> 


1/3 4 


2 1 


2a dt 


for the filter in 


y 


where t = nT = t/(2a) . 


The function is plotted in Figure II— 2, From the figure it appears that a length 
2aNT of 8 to 10 should be adequate, or for 10, 

NT = (x 2 = 10) /2a = 5/a . 


Figures II-3 and II-4 illustrate the amplitude and phase response of three delay lines 
T » .1707/a, N - 32; T = .09/a, N = 64 and X = .09/a, N = 128. The magnitude of 
the error is shown in Figure II-5. The number of taps shown are integer powers of two 
for convenience in the fast Fourier transform. The rule of thumb filter, 

T = .17/a, N = 32 appears quite reasonable to an octave beyond the corner, after 
which the phase response falls off quickly. The second filter, T = ,09/a, N = 64 
is approximately the same length as the first, but with twice the resolution. 

This filter holds the correct phase response to two octaves beyond the comer. 

The third filter, T = .09/a, N - 128 has the same resolution as the second but 
has twice the length. It shows very little improvement over the second filter. 
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5. Conclus Ion - A model for turbulence effects on propagation has been reviewed, 
and a digital representation designed. The design is summarized in Figure II-6. 
The magnitude of the amplitude error is less than 1.5% over a decade of frequency 
beyond the corner. 

Reference 

A2-1 - V. I. Tatarski, "The Effects of Turbulent Atmosphere on Wave Propagation", 
NSF TT-68-50464, 1971. 

A2-2 - C. Abramowitz and I. A. Stegun, "Handbook of Mathematical Functions", 

NBS App. Math Se. 55, Dec. 1965. 

A2-3 - J. D. Markel, "Z Transform Applications Using Digital Computers", 
Engineering Technology, Dec. 1968, 




FIGURE II-1 
IMPULSE RESPONSE 
OF (a/(s+a)) 4/3 
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FIGURE I 1-2 

LINE LENGTH APPROXIMATION 


AMPLITUDE RESPONSE 
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Th(nT) = Yr _1 (4/3)(nY) 1/3 exp (-n Y ) 
Y = .09 
g = .9807 

FIGURE 1 1 -6 

TAPPED DELAY LINE FILTER 



APPENDIX III 

QUANTIZING STRATEGY FOR INPUTS TO VA DECODERS 


1. A Viterbi Algorithm (VA) decoder can function with either hard or soft 
decisions on the input symbols. Typically* the use of 8 level soft decisions 
on the symbols permits the decoder to operate satisfactorily with ^2dB less 
Efa/Ng than if hard decisions are employed. The normally used quantizing 

strategy results in quantization level spacings of about 0.5 to 0.7 times the 
standard deviation of the channel noise. This memo explores the derivation of 
a rationale for choosing the spacings based on the mean and variance of The 
received symbols. The selected spacings are found to vary between ^0.58 and 

0. 64 for means between 1. and 2. sigmas. (3 to 9dB Ej j /N q ). These results are 

equally applicable to a sequential decoder. 

2. Theory of VA decoder operation 

The VA decoder is a device which chooses the most likely state sequence 
of the encoded data from the set of possible state sequences. The formal 
statements of this principle are adequately covered in the published litera- 
ture. The basic operating technique is to associate a path metric with each 
state sequence which allows the decoder to select the maximum likelyhood 
state sequence, based^ on a set of observations. The path metric, is formally 
defined as -In Pr{X|Z}, where X is the set of observations, and Z is the state 
sequence. For a memorly less, binary symetric channel, with Gaussian noise 
statistics, the path metrics can be simply the sum of the observations, weighted 
by the expected observation, given the state sequence for the specific path 
metric. In simpler terms, if M(= +1) is the expected symbol at a given point 
in time, then ^(i) = ^(1-1 , where ^(i) is the Kth path metric at time 

1. The state sequence with the smallest path metric is the maximum likely- 
hood state sequence. 

Any two state sequences will have a number of expected symbols which have 
the same sense, and a number of expected symbols with the opposite sense. 

If one of these state sequences is the correct state sequence, then D, ^(i)= v j < (i)- x ^(i) 
will be less than zero unless ^iX-M^cO, where the set, i, includes ,J 

only those observations where the state sequences differ in interpretation. 

For the correct path, E{X.M^} = n x » 

Clearly, if X is a normal random variable, then^A^ j(i) is also a normal 
random variable, with a mean n A = Nn , and o. = No , where N is the number of 

observations contributing to the decision. Then, at any point, an error is 
committed only if A^ j(i)>0, which would occur with probability QCn^vTT/a^) 

given N pertinent observations. While this cannot be used to directly predict 
the decoder performance, it can be used to predict the degradation in perfor- 
mance when X is not normally distributed. 
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3. Quantization 

When the symbols are quantized prior to decoding, the X's are no longer 
normally distributed, since they take only integer values, and finite range, 
linear quantization results in significant population groups in only a few 
values of X. Consider the quantization strategy shown in Figure II 1-1. We 
have arbitrarily assigned integer values to Z, the quantized estimate of X. 
If X is a normal random variable with mean n and unit variance, then the 
probability distribution of Z is given by, x 

Pr{Z = -7} = Q(n +3a) 

Pr{Z = -5} = Q(n^+2a)-Q(n x +3a) 

Pr{Z = -3} = Q(n +a)-Q(n +2a) 

A. A 

Pr{Z = +5} = Q(n x -3a)-Q(n x -2a) 

Pr{Z = +7} = 1-QCn -3a) 

X 

p p 

where Q(X) = 1 /" exp-{t-n ) /2a dt and "a" is the quantization level. 

a 72V x XX 

X 

It is rather awkward to compute the probability that A>0 directly, since 
the probability distribution of a is the result of N convolutions. A more 
useful approach is to use a Chernov bounds. In this case, we derive these 
bounds using the Laplace transform of the probability distribution of Z, pro 
ceeding as follows: 

♦ 7 (S) = E 8 Pr{Z = 2J-9}e S ( 2J “ 9 ) 

L J=1 


Since convolution in normal space is simply a product in transform space 
the transform of the probability distribution of a is, 

VS) - [* Z (S)] N 

Let v A (S) = In $ A (S) 

¥ a(S) = Nln[ I 8 Pr{Z = 2J-9}e S ^ 2J-9 ^] 

J=1 

N z 8 (2J-9)Pr{Z=2J-9}e S(2J " 9 ^ 

J=1 

E 8 Pr{Z=2J-9}e S ( 2J “ 9 ) 

J=1 

Now, the Chernov bound can be expressed, 

PrU>o} < e T4 < s >| 


* A (S) = dT A (S) = 
~dS 


%'(S)=0 
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Thus, we need to solve for the value of S such that, 

E 8 (2J-9)Pr{Z=2J-9}e S t 20 " 9 ) = 0 
J=1 


This can be solved numerically, and a value for Pr{A>o} found for any 
however, that 

v,(S) N 


N. Note, 


v-' / ■ 

Pr{a>ol<£e 3, 

Thus, we need, only solve for the case H-l. Jgd «■£ the of 

x’ 


minimized Pr{A>o} For the case ^ 

interest the variation of Pr{A>o} with "a" is shown in Figure IH-2. We can clearly 
see that’varying "a" over rather broad limits does not change the probability 
value greatly, however, the optimum is clearly in the vicinity of a - - 62< Y 

Figure 1 1 1-3 illustrates the same information over a wider range of mean to 
standard deviation. 


Now we note that there is a hazard in using any bound to draw c ^ c ]^sion 
of this’nature, since varying the parameter H a" may result in v ®|' iatl0l ? s .]” ^ e 
relative tightness of the bound which could be on the order of the variations 
in the results Thus, if our conclusions were at significant variance wit 
emperical data’, it Sill be necessary to use a more sophisticated approach to 

eliminate the uncertainty resulting from the i us %S f t £ U !!lL 1c ?TSsilts and* 
that- thp pmnprical data is in close agreement with the predicted results, ana 
that the optimum is rather broad, suggests that the increased sophistication is 

unnecessary. 


4 ‘ It n has S been shown that the criteria of minimizing the Chernov bounding value 
of PrU>o) leads toa choice of quantization level of n,0.6o x . which is halfway 


between the emperical ly determined values of 0.5 to 


1 1 1— 3 




FIGURE III-l 
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FIGURE III-3 

EFFECT OF VARYING MEAN ON 
CHOICE OF QUANTIZATION LEVEL 
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APPENDIX IV 
LIST OF SYMBOLS 
(Input Symbols) 


ASOFT 

Soft decision quantization 

reference 

BIF 

IF Bandwidth 


BL2 

2 B lq (AFC) 


DELF 

Delta frequency 


DYNR 

Dynamic range of AGC below 

"minimum" signal 

ENODB 

Energy to noise density in 

decibels 

FLO 

LO detune 


FQDPRT 

Frequency Doppler rate 


FRQDOP 

Frequency Doppler 


11 

Random generator number 


111 

Random generator number 


12 

Random generator number 


122 

Random generator number 


NSLICE 

Phase of bit synchronizer sample data filter 

SDAMP 

Gam of amplitude scintillation 

SDPHA 

Gam of phase scintillation 


SVRHZ 

Sweep slope in Hertz 


SWNM 

Sweep minimum 


SWMX 

Sweep maximum 


TL 

Threshold 




(All Symbols) 


A Signal voltage 

AAMP Normalized amplitude scintillation 

ACOR Scintillation normalization 

ACQ Number of true acquisitions 

ACT HI A cos (TH1 ) 

AFAL Number of false acquisitions 

ADB A dropped bit 

ADOA Average number of bits in a dropout 

ADOS Standard deviation of number of bits per dropout 

AGC AGC gain function 

AK AFC Loop gain 

AKOU KDATA at acquisition 

AKBMIN Bit synchronizer AK minimum 

AMIN Minimum A to AFC 

AMIU2 AMIN/SQRT (2 ) 

AMIN3 AMIN for dynamic range 

AMISS Number of missed acquisitions 

APHA Normalized phase scintillation 

ASOFT Soft decision quantization reference 

ASOFTI Soft decision quantization slope 

ASTH1 As in (TH1 ) 

AVEES Average ES 

A5 Baseline connection filter gain 

BAMP Amplitude scintillation approximation intercept 

BIF IF Bandwidth 

BL2 2 B lQ (AFC) 

BL2B 2 B^q (bit synchronizer) 

BPER Bit period 

BPHA Phase scintillation approximation intercept 

BRATE Bit rate 

CAMP Complex amplitude into BIF 

CBINI Hits in SBINS1 
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CBIN2 Hits in SBINS2 

CBIN3 Hits in SBINS3 

CBIN4 Hits in SBINS4 

CB1 Bit synchronizer Z transform constant 1 

CB2 Bit synchronizer Z transform constant 2 

CB3 Bit synchronizer Z transform constant 3 

CCL Control function for scintillation subroutines 

CCLT Control function for TRIG subroutine 

CCCUT Output of last IF amplifier 

CDTHX Cos (DTHX) 

CDTH1 Cos (DTH1 ) 

CFILT IF filter output 

COUT CCOUT 

CSCIN Cos (YP) 

CTHX Cos (THX) 

CTH1 Cos (TH1) 

CTH2 Cos (TH2) 

Cl AFC filter Z Transform constant 1 

C2 AFC filter Z transform constant 2 

C3 AFC filter Z transform constant 3 

C4 Normalized AFC VCO gain 

C5 AGC Z transform constant 1 

C5DP Doppler Z transform constant 1 
C5MH Baseline Z transform constant 1 

C5MS Sampling Z transform constant 1 (word "sampling" is synonomous with 

C6 AGC Z transform constant 2 word "slice" in computer listings) 

C6DP Doppler Z transform constant 2 

C6MH Baseline Z transform constant 2 

C6MS Sampling Z transform constant 2 

DELF Delta frequency 

DES DT*SVR 

DT Delta time 

DTHX Delta THX 

UTH1 Delta TH1 
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DYNR 

E 

EDOP 

EDOPR 

EDOPRX 

ENO 

ENODB 

ENOM 

ENOMIN 

ERR 

ERRM 

ES 

EV 

EVMAX 

EVMIN 

EVO 

EV1 

EV1C 

EVIO 

EVIOC 

EO 

FBAMP 

FBPHA 

FDIFF 

FDIFFB 

FIDO 

FLO 

FQDPRT 

FRQDOP 

FVCOB 

GAM43 

GAM43A 

GAM43P 

G1 


Dynamic range of AGC below "minimum" signal 

Voltage into AFC filter 

Doppler voltage 

Doppler rate voltage 

Doppler plus Doppler rate voltage 

Energy to noise density 

ENO in decibels 

Minimum ENO 

ENOM in decibels 

Error count 

Error spacing matrix 

Sweep voltage 

Voltage from AFC filter 

Maximum sweep voltage 

Minimum sweep voltage 

Last EV 

Voltage after Doppler addition 
Voltage into bit synchronizer 
Last EV1 
Last EV1C 
Last E 

Amplitude scintillation filter corner 
Phase scintillation filter corner 
AFC pull in 

Bit synchronizer pull in 
Float (IDO) 

LO detune 

Frequency Doppler rate 

Frequency Doppler 

Frequency VCO bit synchronizer 

Scintillation filter tap spacing target 

GAM43 for amplitude 

GAM43 for phase 

First IF amplifier gain 
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G2 

Second IF amplifier gain 

I 

Dummy index 

I AMP 

Amplitude scintillation counter 

IAMPC 

Phase scintillation counter 

I B IT 

Phase of bit synchronizer 

I DATA 

Current bit 

IDATAO 

Last IDATA 

IDATAW 

Last IDATAO 

IDLBIT 

Number of output bits before error count 

IDO 

Number of dropouts 

I DON 

IDO counter 

IERRSP 

Error space 

IH 

Histogram bin 

IHOLD 

In phase bit synchronizer level 

IHOLD1 

Last IHOLD 

ILOC 

Counter, bits to lock 

ILOCF. 

Final ILOC 

IMS 

Hard decision from sampling filter 

INIT 

Control function for BIFILT subroutine 

IPHA 

Phase scintillation counter 

IPHAC 

Amplitude scintillation counter 

IX 

Bin of soft decision 

IXM 

Last IX 

11 

Random generator number 

ns 

11 to start 

in 

11 for scintillation 

ms 

111 to start 

12 

Random generator number 

I2S 

12 to start 

122 

12 for scintillation 

I22S 

122 to start 

J 

Counter for samples per bit 

K 

AFC VCO gain 

KB 

Bit synchronizer VCO gain 



KDATA Counter of IDATA 

KOUNT Counter of IMS 

NACQ Counter of NACQT trials 

NACQT Number of acquisition 

NBIT Phase of IBIT to dump in phase integrator 

NBIT2 Phase of IBIT to dump quadraphase integrator 

NC Quadrature noise 

NRUNS Number of bits per run 

NS Quadrature noise 

NSLICE Phase of IBIT to sample data filter 

NSPB Number of samples per bit 

NO Noise density 

N43 Number of scintillation filter taps 

OMC Baseband equivalent of IF in radians 

PACQ Probability of acquisition 

PFAL Probability of false acquisition 

PE Probability of error 

PEAK Peak voltage from sampling filter 

PHIB Steady state bit synchronizer phase error 

PHIDEG Steady state AFC phase error 

PHISS PHIB in radians 

PHISSB PHIDEG in radians 

PI Pi 

PI2 2*Pi 

PMISS Probability of missed acquisition 

PKAV Average of IBIT when PEAK occurs 
PKSD Standard deviation of IBIT when PEAK occurs 
QERR Voltage to bit synchronizer filter 
QERRO Last QERR 

QV Voltage from bit synchronizer filter 

SBINS1 Probability matrix for -1-1 input 

SBINS2 Probability matrix for +1-1 input 

SBINS3 Probability matrix for -1+1 input 

SBINS4 Probability matrix for +1+1 input 
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SDAMP Gain of amplitude scintillation 
SDPHA Gain of phase scintillation 
SDTHX Sin (DTHX) 

SDTHO Sin (DTHO) 

SDTH1 Sin (DTH1 ) 

SIG Normalized noise 

SKOU Standard deviation of acquisition at KDATA 
SLOA Amplitude scintillation approximation slope 

SLOP Phase scintillation approximation slope 

SLSIG Standard deviation of sampling output 

SLSUM Mean of sampling output 

SSCIN Sin (YP) 

STDEV Standard deviation of measurement 

STHX Sin (THX) 

STH1 Sin (TH1 ) 

STH2 Sin (TH2) 

SVR Sweep slope 

SVRHZ SVR in hertz 

SWMN Sweep minimum 

SWMX Sweep maximum 

TAUB1 Bit synchronizer filter constant 1 

TAUB2 Bit synchronizer filter constant 2 

TAU1 AFC filter constant 1 

TAU2 AFC filter constant 2 

TAU3 AGC filter constant 

TAU4 Sample filter constant 

TAU5 Baseline filter constant 

TAU6 Doppler filter constant 

TEMP Temporary duirany 

TH2 Phase angle from AFC VCO 

TL Threshold 

TSYNC Period of bit synchronizer VCO 

U Voltage into AGC filter 

UMH Voltage into baseline filter 
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V 

Voltage from AGC filter 

VO 

Limited version of V 

VDP 

Voltage from doppler filter 

VHIST 

V histogram 

VMH 

Voltage from baseline filter 

VMS 

Voltage from sampling filter 

VMSS 

VMS at sampling time 

VMSO 

Last VMS 

VSQRT 

Square root of VD 

WN 

AFC loop natural frequency 

UNB 

Bit synchronizer natural frequency 

X 

Dummy 

XI 

In phase integrator voltage 

X I B IT 

Random starter for IB IT 

XIO 

Last XI 

XIOO 

Last XIO 

XISIG 

Standard deviation of in phase integrator 

XISUM 

Mean of in phase integrator 

XQ 

Quadraphase integrator voltage 

XQO 

Last XQ 

YA 

Scintillation amplitude 

YP 

Scintillation phase 

ZETA 

AFC damping factor 

ZETAB 

Bit synchronizer damping factor 
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APPENDIX V 
ERROR RATE 
COMPUTER SOFTWARE 


V-l 



a ooao3- 

1)110003 

0 00003 - 

300003 

000003- 

000003 

000003- 

000003 

000003- 

000003 

OQ0Q 8-3- 

000003 

00000 3- 

000003 

0 00 00- 3 - 

000003 

i n n n 3 

U u V U U 

000003 

000003 

000003 


0 000 03- 
00000? 
0 00 0 1 0- 
000013 
00001 4 
00001 ? 
0 0 00 go 
090022 


000022 
00 00 25 - 
000026 
000027 
000030 
000031 
000032 

I'.inuri 

"v VU U vO 

000034 

00 00 35 

000036 

a nn m-7 

V UU U vr 

000041 
0800 42 - 
300043 
0 000 -1 4 
0 00045 
0000 4 6 
900047 
000051 
000052 
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PROGRAM C ARL( INPUT, OUTPUT, TAPE5*INPUT,TAPE6*0UTPUT> 
COMPLEX - CAMPT 4 ) ,CFILT ( 4 ) , COUT , CCOUT 


REAL NS,NC,NQ,K,KB 

Dl f ^ 'STPM S9TMS1 ,, P n tl i S 3 1 % T ft» T S«TMSh < & r 

01 MENS I ON AOB (IOO) 

— DIMENSION F1(8)>F2W 

DIMENSION £RRM(100) 

DATA FBAMP , FBPH ft ,SP A MP , SQPHA , NQ/2 ., .2, . 23,. 4 7,1./ 

DATA GAM43,N43/.09,64/ 

D ATA EN0D8, BRAT£,EN fl MIM/ 6 . , 6 8 » , 9 . /— ; 

DATA OELF,NSP9,NRUNS,BIF/62.*40,1Q, 1000./ 

D A TA gLZ y PHlBEG^OIFF/lT B.i lO . *3 4 2 . ^3 / 

DATA 8L29,PHIB, F0IFF9/1 ., 10 • * 5 •/ 

DAT A I - l -,- I2 , 111 »- T22 * NSLICE - /1 * 1 , 1*1 , 2 4 / 

DATA FRQOOP*FQOPRT,TL , AGC/0 . , 0 . ,-2. 6, 1./ 

DATA TAU3 ,ASOFT/ . 07957 , . 2 4 7 4 / 

DATA TAU4*TAU5,TAU6 ,A5, FLO/. 00180 8, .01808 , . 10603* .5,0./ 

DATA SVRHZ , SWMX , S H HN/ 2 . 2E 3, 28 . 15E3* 6 . 8 5E 3/ 

DATA IDL8IT/20/ 

DAT A OY N R/O « ■/ 



NAMELIST/INI/ NRUNS, NSPB*EN0 09»0ELF * BIF,EMOHIN 

» * BL2 > PHIOEG , FDIFF , BL2B , P H I B , EDIFF B 

+ » FBAMP , F9PH A, SOAMP, SDPHA,GAM43»N43*NSLICE 

» * FLO , FRQPOP * FQPPRT * ASOFT* II , 12 * 111 , 122 

+,SV*HZ,SHMX,SWMN,TL,AGC 

» * I B IT * TSV NCfIDLBiT* OYNR 


’f,TAU^ 

PI^.»ATAHC1 . » 

PI2=2.*PI 

ZETA-SQRTt0 . 5) 


ZET AB=ZET A 


CCL T*fl » — 

CALL TRIG (CCLT ,PI *SSC IN *CSCIN ) 


C CLT -jn 

CONTINUE 
RE A D INPUTS - 
READ(5,IN1> 

NSLICE-2 6 

CBIN1=Q. 


0 8 IN2~0 -r- 
CBIN3=G. 
CBIN 4 -0. 
ILOC= 0 
ILOCF— 1 
100=0 
IDOM--1 
FLO=10. 

P YNP=3 . 
11 = 1 

12=1 


111=1 


122-1 

FQDPRT=19 .8 
NRU N S=2Q 0£— 
TL=-b .5 


S VRHZ=d - ^ 
AS0FT=.3 
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oooos3 - 

a -JO 054 
430356 . 

00006a 

00 03&1 - 

J 00 064 

0300o6 

300367 

000071 

006073 

3001375 

00a076 

303110 

JuOiOi* 
003111 
000117 
00 a 1 21 
306126 
003126 
0 0 0 1 30 
0 00 1 35 
000167-.- 
000141 
0 fl« 1 52 

000157 
300162 
0 00 167... 
a a 0 1 7 1 
000176 
0 0320a 
000202 
00020- 
6 0 0 2 v ^ - — 

3 30 211 
000216 
3 0*i215 

U GO 2 2a 

3 3 3 2 25 
j 002 27 
3 J C 2 34 
000240 

3 0 a 242 

J 0 0 24 c 

0 00246- -- 

0 j a 2»1 

000253 
303255 
0 u 02-66 ■ 

a J 0 26 1 
0 00 266 
300265 
•j 0 0 267 

a002 7 2 

u 0a27^ 

033277 


- - I-OL-BIT=lQG 

AS0FTI=2«/AS0FT 
- 00 199 1=1,6 
00 19 8 J=l, 8 



S3IMS2(T»J)=0 • 

S3 1 NS 3 (I, J) =0 , 

S3INS4 Cl , J) =0 . 

198 CONTINUE 

199 CONTINUE 

iXM=0 - - -- 

CALI RAN0A(I1,I2,XI8IT) 

IRIT=INT(39.*XIBIT+1.) 

C CONVERT INPUTS ( ENu , 3PE R , PHIS S , PHISS 9) 

tNu=10*** (ENd03/10. ) 

E NO M= 1 0 . * * C EN J 'll N/ 1 0 . 3 

. . 3R£E= 1 • / ORATE 

3 T = e°€R/MSPP 

PHISS = PHI0FG8PI/18o ♦ 

Q HlSSS=PHTB*PI/ieO. 

AMIN=SQRT (FNaM*NO*BRATE> 

AMIM2 = AMTN/ 1. 4142 

- S6C = S CRT (AHIM-2) • - ■ 

AMIN3=AMIN2/SQR T C1C .♦ * ( DYNR/1 0 . ) ) 

. ASSORT <EN«*N* ♦ORATE) 

0 , CALCULATE LOOP PARAMETERS ( T A Ul ,T AU 2 , K , T A JB 1 , T AU9 2, K9) 

WW=SL 2/ (2ET A+0.25/ZETA) 

AK=PI 2* F 01 PR/ SI NCFHISS3 

T AUl r AlC/ {WN*WN> - - 

TAU2=2.* ( ZETA- j ,5/ ( WN*T AU1) ) / WN 

K-AK/AN-IN 

E00F = PIZ* FRQOOP/K 
E00F p = PT2 # P0DPRT* DT /K 
cDOPPX=EOOF-*-Ennpp; 

- S-V1LA X- RJ t* S W« X-/ x 

EVMIM=-=I2»SWMN/K 

SVP=PI2*SVRH7/K 

MM3 = 9L23/ (ZET A8+C.25/ZETA9) 

AK = PI2TFr)IFFP/SIN(PHISSB) 

T AUB1 = AK/ (HNB*NNB> 

XAlJ32=2^*-CZET AO-O .5VC WNB^TAURl) 3/KNR 

AKQMT N=PT 2* (0 EL F/ 2 . ) / ( K*1 . 4 12 5 ) 

KB=AK/AKOMTN 
C OlFINE CONSTANTS 

Cl=l. -QT/TAU1 
C2 = T A U2/T AU1 

.. 03=01*02- 1» — • 

Cf- = nT*< 

C 5 =■ 0 T / TA U 3 
C6=l. -C5 
-jrHl=PI^OELF»DT 
SOT HO = S IN (DTH1) 

. . ^OT-Hl-COS COTH1) 

OTHX=PI2*FLO fc OT 
SDTHX*SIM (O t HX) 

COT HY = C OS ( OTH X) 

Cal=l ,-BPER/T AU91 
C B2 = T AU92/T AU91 


V-3 



0 Ou 3u 3 

nn m f) c 

— ww v w y / 


82-1. 


C5MS= CT/TAU4 
-£6-NS=-W=£-5MS 


Q0Q.3Q6 

-00*340- 

oou3u 

-04 0313 


C5HH=DT/TftU5 
CoMH=-l-.»r-5W4 • 


0 Ou 31% 


000321 
■3 -& 03 22 - 
000323 
0 0 0 32V- 


C50P=0T/TAU6 

(^-QR=l-.--&5£P- 

SIG=a .5*$QPT(NO/OT) 

I44.I4I-A4-IZE --tOGP-S- < T H2 * £ 8 **I »XQ t GV » QER-RO-r£-RRi- 

SLSUM=0. 

St * 


000325 

-aa&i30- 


XISUM=0. 

XlSi<i=0<. - 

N8IT=9PER/9T*.5 

US£X2^U&IT/-2- 

T SYNC- 8PEP 


300335 



CftMPf 11 a 

CAMP(2) = (0.,2 .) 

000342 


v y » t — 

CFItT ( 1 ) = (0 » » 0 * > 

-0-30 344 

000347 


—ii t i-L 1 \ ci — 'U • ♦ -u ■• f — 

CFItT <3)= <0 ., 0 .) 

-0-0-035-1 — - 
000352 


loi A “ 

0MC=?I*9IF 

ID ATA ~-l - — 

000355 

ft ft 4 

PRINT 20 1 » ENii DB»DELF»9IFjBL2» NSLICE » ASOFT 

crtouATH y * pki n n n- * . Fk * 1 «X 6 » !♦* * F 6 * Q + — ■ 

il 13 U J f ? CM^L r 7’^ - » * ~ * ' 

+*■ » BL Z-* » F 6 • 0 t* » NSLICE=*»I 6 , * , ASOFT=* i .4) 

3 BTHT ofto oiauT .cunv. cumm.fi n.cpnnnp.FnnPAl 

000415 

2C2 

FORMA T (IX »* SVRHZ=* » £12 • 4»* » SWNX=* » E12 . 4 » * i SWMN=*, E12«4» 
j.* . tri n-*.Li 2 . u-,*' , rgQnoPjj+>£12*4 r » , FQ0PRT=» . £12 ,4) 

000415 

ft ft ^ 

PRINT 2 03 ,TL,DYNR,SOAHP,SOPHA,Il f I2,Ill»I22 

eft n u a t M v *T 1 Ng=» * rgi-1 * 

y il n 1+14-1. CMM r u n riw 4 - \ o> ■» 1 — *-4 t ’ 1 ' 

+* »SDAMP=**F 6 .3» ¥ t SOPH A=*»F 6 .3»* » 11 > I2» 111* 122* ,416) 

, « « n o T kit 9 HL. UM. I^.TAIM T &11?^UMQ - K-fl - T All fll -T A — 

M UU44J. 

000465 

204 

F 0 PMAT< 1 X,*WN-*, F 6 .lt* ,K=*,F6.1,* ,T AUis*,E12.4f * ,TAU2=*,E12.4 
.uM Be *.K.H.4 -Kp=» t £J»^ r *--^I-AUai=y f E12.4i7 .TAU82=*tE12.4) 

000465 


I ERRS P-0 

Q4Q4b6 

00047? 

1 ^),, 

AOB (I > = 0 • 

000474 

- ■ C CM- 

!j = 0 . 

z u f j — n . — 

0 QQ47 5 

000476 


b-V o - u ♦ — 

TH2 = 0 . 

PTU9- < — 

000500 


nc-lr+ 

STH 2 =a. 

aoofcai — 
000502 


STHX=0 • 

ci *i — r 

OQpy 3 — 
000504 

n n ri c nc 


EV= B» 

VT* 11- 

y u Op Up — 
000506 

ft nncnT 


Al^U ♦ — — — 

XQ=C . 

000510 


QEPR0=Q. 

n - 4 ¥ —— — 

~yy ypil — 
000512 
0 005-13 


b 

G 2=1. 

4*^1 — 


000514 


E$=0. 




JQu5i5 sia-isa* . 

J 005 16 CTHl=i. 

0045.47 JJ£SaOTAS\|R - 

000521 ERP-G. 

4455-22 -EgROHPsQ... 

0 0 05 23 K OUNT = -l 

-0 04524 IHOCO.S4 . ..... 

000 e i25 T H 0 1 n 1= 0 

• 3 00524 4MS=£4/-EQ0P . 

000530 EV1 0= VMS 

3 305-31. -ElLLOCsVJdS.. - .. 

0 0,05 32 VHH=G. 

304533 . UMH-sO-. . .... 

000534 VHPrO. 

334535 -XlOOsfi. .. 

000536 X T 0=g . 

3 3.45 37 J?. NSP 3 .. .. . 

0J054C °KAV=0. 

0045 41..— PK$ D=. 0 • 

r TUP PUL ENCE IMITIALI ZATION 

000542 ACOR- SQRT 12 7 5 “! •) _ .. 

' 000551 AAMP=PI2*F BAMF/ACOR 

4.04553- A£HA*£I2*FaeH A/AC OR. „ 

0 0 0555 I AMP = I NT ( ,E«-GAH4 3/ ( AA MP*0 T) ) 

44u561 IPRA^INTU5+GAM43/{ APHA*OT) ) 

0 0056a G AM43 A=OT*FL0AT (IAMP)*AAM° 

000567 GAMA3P-. r )T4FLO-AT.(IPH.AJ *APHA 

000572 CPL=0. 

040 5X3. — C AU.-. 450 5JUJJ4 3 , G AW* 3A » S DA M CC L fZ AHP t SL O A* a A M P , 1 1 1 , 1 22 ) 

3 0 u 6 u 4 CALI BL03F<N*3,GAM43P,SOPHA,CCL,IPHA,SIOP,3PHA,I11,I22> 

4 3 j 61 5 - €0.-1. 

030617 TAMPC=0 

043620 — I PH AC *3 

C SEL FC T DATA (IDAta) 

44-0521 tCD A T A =4- - _... 

000622 a CON T T NtJE 

si 036 22 - ... ..IAMF.CsJ.AMFC *1 

300624 IPHAC = IFMAC+1 

044625 — ... I -Ft IA.MPC* GT . I AMPJ GO TO 10 

00J63C 13 IF(I°HAC .GT.IPHA) GO TO 11 

0345 3*. GC-TO. .12 .. .... 

303634 IQ IAMpr = l 

040635 CALL BL OQA (N* 3 , GAM4 3A ,SOAMP »C CL ,IA HP, SI OA , BAMP, T 11, I 22 ) 

3 0 0 6+6 GO TO 

000647 11 ... I PH AC - 1 

00365C CALL BL 03P ( N4 3 , G A M4 3P , SOP hA , C CL , IPH A, SL O 3 , 0PH A , 111 , 122 > 

.y 4 0661 12 CChtXI-61U£. . . ..... 

000661 YA = SIOA*FLOAT (IAMPO + 5AMP 

444665 XE=SLOP*FLOATlIF>HAC)«-ePHA 

340670 CALL TRIG f CCL T , YP , SSC IN, CSCIN) 

440673 . ... XF-d^ir^HSPPlCO TO 1 - ....... 

300676 J-3 

445676 XC.AUUlslHA.T3U . — 

04J700 IDATA4=IDATA 

000701 ... CALL- -OAT AS (1-1 AT A) ...... 

s,<)Q7u2 SDTH1=POTHO*IDATA 

440705 IPlV^GTVr-LJGn. - TO 1 

30u71C ^KAV=PKAV+ d EAK 


V-5 



000712 

003714 

30-714 

r 

00 — 716 1 

3G-7lb 
0007 ?2 

3 00*24 
0007 26 - 
030731 

3007 03 - 

300*34 

0(107^*4 

000*34 

-3076b 

000760 

000763. 

000765 

0GC767 

r 

000770 

00 10 j? 

00101-2 - 

J01-2E 

301-22 

r 

041026 
001027 
001033- 
001036 
0 0 104m 
0-10 ♦ 1 

c 

- U 10 46 

OOiQ^oi 

001 - 61 

001367 
301066 
001072 
301074 
0-1077. 
301102 
00 11-* 
0-1137 
001112 
001120 
00112-4.- . .. . 

- - 1 1 27 
001130 
001132 
001133 

0-11 34 
0 011.36 
0011.2 

G 

0-1144 

031146 

OJll.* 


PKSJ>PKSO + PEAK*P£AK - - - . - 

PEAK = Q * 

_ JCaATAslCOATA*!-- - • - - - 

DEFINE INPUTS (TH1,NS,NC) 

.CONTINUE — 

TEMF = STHY*COTHX<fCTHX*SDTHX 

. CTHX=CTHX»CDTHX^STH3(*SDTHX . — . - 

S7HX=T£MP 

. TEHF=STH1*CDTH1 + CT^1*S0T6<1 - - 

ETH1=CTH1*C0TH1-STH1* sdthi 

. -,-SXW 1 - T-E-MR. — 

4STHl^l.4l42*A*tl.*YA)*(STHl*CSCIN4CTHl*SSClN) 

ACTH1 = 1.-4142*A* (1 „**A.) *.(CTH1*CSCI N^S7*4i*SSCIN4 

CALL NRANQH (11,12, X) 

NC-X* STG ... _ . 

CALL NRANDM (II, 12 ,X ) 

MSsX* $14 •— ..... 

I 3IT=TPI T +1 

J=J+i _. 

LOOP POUATIOMS 

CANF (A)=CMPLX (ASTH1 +NC, -ACTH1 -NS) *GMPLX (CTH2^ -ST4424*G1 
CALL 01FTLT (CANP,CFILT,CNC,DT,INIT) 

CCOU-T =G2*GFTL T-(-4 ) 

COU T =CCOUT 

V=C6*V+C6*U 

ACC 

VO = V . -- 

IF (ADS(VD) . LT.AMIN3)V0= AMIN3 

- XF-( V0-, NF . G. 4 C AELr- -AS-QftTF 1 -O-Ot VS QRT 4 

Gi- agc/vsort 

r,z~r-i - . _ 

U- AINAG(COUT) *CTHX-REAL (COUT) *STHX 

AFC ... . 

E= PEAL (COUT)^CTHXfAIMAG(COUT) *STHX+ES 

.._c.y=ci.»Ev*€2*£-G3*E0 - - - - 

c V1=E V 0 “EOOPPX 

IF<tV.GT.tVMAX»OES=-OES - 

IF (EV.LT.EVMTN) OES=-OES 
EDOFR X=EDOP=>X+£DOPP 
T H2 = T H2+^^ + £Vl 

. -T-H2— AM OG ( T 2) - - 

CALL T°IC (CCL T,TH2, STH2 ,CTH2) 

VNSO=VMS .. 

VNS=C6MS* VM S+ C5 MS # E VI OC 

IF(ABS(OMS) .GT. A9S( VMSO))P£AKsFLOAT (IDIT) ...... 

OOP =C6DP* VDP+C5 0P*E VI 0 

Vi-C --E Vl-VHH- V4D- -■ • — - — 

XI-XI+EV1P 

XQ= XO + EV1C . . ... 

E V J = F V 

£■-• = £ . . .. 

EV10= E Vi 

— IF IV T-&T-.-7 L 4 EF^E-S + BE S - 

Evmc =evic 

PIT SYNC -TIMING — 

IF (I3IT .NE. NS LICE) GOT 081 

IH-S--1 

TF( VMS.GT.O .> IM^=1 


V-6 





-tfM SS - VM S-- 


001152 81 CONTINUE 

001154 XQG=XQ/FLOAT INBIT) 


2 


OU 11 57 

n a a a c.'y 

2 

■ 

I FCTBIT.LT. NBIT) GO TO 4 

tc /i/rutMT it T m DTT1 (tnTnO'i .. - 

001161* 

*n ,i a 4 ay 


lr UkUU il-f- *51 t W+UTi" - — 

XISUM=XISUM*A 9SUI) 

’ U U i-t-Or 

001171 

n n 4 4 94 

91 

CONTINUE 

rum in -rum n — - - 

001173 

flrt4 4 7(i 


IHOIO= -1 

u y-i i f 4 — 

001176 
n ii 4 onn 


XI00=XI0 

UU lcUu — 

401202 
nn 4 


VMH=(C6MH*VMH+C5MH*UNH) *A5 

1 1 MU— 11 . _ 

001206 
n fl 4 94 9 


IFCIH0LD.NE.IH0LD1) UMH=XI00*Xl0 

% J-Q A - . 

UUliifc 

001213 

fliH 04 ? 


QERR-XQO* ((IHOLO-IHOLOD/2) 

mi— ra< »nwif'i)')«nrDD. ro 3*flFRS 0 

4/ li l iHrt' 

001224 

nn •! 


QERRQ=QERR 

TRTT-fl - 

001226 
n n 4 


FVCOB=BRATE*KB*OV 

T g VMr-TSVWt -N BTT*PT m/ Ftf COB 

UUIc 

001236 
n n 4 ■»/. 4 


NBIT=TSY NC/9T+.5 
MBIT2-NBTT/2 

001243 

n n -i 


IF (V» GT.TL) GO TO 95 

ILQCF-ItiQC 

001247 

n n i 9 #. 9 

402 

CONTINUE 

KQU NT = K0UNT#-4 

■IHi-l'CTf 

001251 

n r» ^ 


IF (KOUNT.LT .I0L8IT) G0TO93 

T c 1 me MC _T HA TAfl ) FRR- FRR-H 

UUi&yi) — 

001257 

r 

1 Xi 1 1 Uj v lie, * 1 yn IRIM t- i>* ““ trvftTi m 

IFdHOLD.NE.IDATAO) ERRO HP-ERRDMP#-1. 

FBROR SPACING - ----- - 

001263 

ft n 4 oc.it 


IFdMS.EQ.IDATAD) GO TO 221 

T tf TCDDCD -1 T.MTrOOep^i 

001270 

n n ■< 


* r 1 IE J ^ ib 1 1 It- l\l\ W « X ■ 

IFdERRSP.GT. 100) IERRSP=100 

c ppMi tcoogpi -EflflM f TFRRSP) tl x 

001276 

ft ft 4 9 7C 


IERRSP=0 

emo? ?? — — 

001277 

naiini 

221 

IERRSP=IERRSP+i 

TONTTMUF 

uvlvvi 

001301 

nn 4 9 ft u 

etc 

slsum=sLsum*abs(vmss> 

U W * V W T 

n 04 int 


OLOlU"9VOiU T f “OiO »n O J 

SOFT DECISION 

T V— TMT 1 *CriCTT» uue- - in 

001312 

11 ft 4 74 C 


IFdX.GT.8)IX*8 

Tp I T y _ 1 T_HTy;1 . 

V tf A<>15 

001320 

11IH JS4- 


1 t 1IA IL ! 1 1 r 1A 

I FdXM »EQ »0) G0T096 

TCiTnfctAu irn aud t da TA A - f T V _ ▼ t t v tvu\ . a 

001334 

n ft 4 ii. c 


in 1UH 1 HR • L U • l«nin/«*UH 1 H IM * .x * w Jill Jx W, n | x Al II * J U1 n^x 1 I A j X ft Tl f P X 

IFdDATAW.EQ.-l. AND.IOATA0.EQ.-1) C9IN1=C3IN1»1. 

rc/Tftmu trn _< tun Tn»tin rn nirnrucoiTV. tiiui-potum / tu Tvui n 

001360 


IFdD AT AH .EQ. -1. AND»I DATA0 • EQ •+1)CBIN2 S C3IN2+1« 

rcitn»t«u rn iii^.TAAIA.A.Cn — SAT-MC XtX V- .tvm» — c*qt upj t tv tvui .4 

001404 

001415 


1 T \ 1 U H 1 M TT • HI* nnu * x u n * n U • Lva < i# iiw w ^ X ft 7 1 a hi g ux 11 l 1 A 7 1 AH| ?1 

IFCIDATAH.EQ.+l.ANO.IDATA0.EQ.-i)C9IN3=CBIN3H. 

IFCI0ATAH>EQ«»1.AND»I DATA0»EQ .♦!) SBINS4 CIX« IXH) ■SBINSb dX,lXH) H 


001427 


IFdOATAW.EQ.+l.fiND ,IDATA0.EQ. + 1)C8IN4=C3INI»*1. 


V-7 





IF(ILOCF.EQ.ILOC) GOTO40 


001444 


801447 


001452 


001455 




031461 


001465 

001 4 65 - 

031466 

n fl i /i 


1471 


001505 
-4415 tt- 
001517 


IffiKVHkXI 


0u 15 36 

001536 

-044542- 

001552 

0 0155 6 

001564 


IF(ILOC.GT.NRUNS) GOTO40 8 

G0T0 9 3 

IF(IDO.EQ.IDON)GOTO401 

rnn-rnftji< 


IF(IDO.GT.10C)GOTO402 


Ihlal.matui 


409 (I DON) = 


k vr ■ ■ ■• a«i>' 


GOT0402 



FORMAT (IX,* ACQUISITION FAILURE*) 


CONTINUE 


ERROR RATE CALCULATION 

OnuTTantr 


PE=ERR/(K0UNT-IDL9IT) 




SLSUM=SLSUM/(KOUNT- ID L 8 IT) 
-S1S T G ~ SQ R T ( SL SI G / (-KOUNT -IQL 
KKK-KCUNT-IOL9IT 


iiipiyg 




FORMAT (IX ,*PE=*,E12 .4,* , STDE V=*, El 2.4* * ,V MEAN=*,E12.4, 

PEI=ERROMP/ (KOUNT-I DL PIT) 

5 PEI- SORT (PEIMl-r -PEI),/ (KQUNT^-EQ LPI T)-) 

XISUM=XISUM/(K0UNT-I0LBIT> 

- X I S 1 & 8 . S QRT4-XI-S46A(X-QUAIT"-IQL-9I T ) » X I S U M*X45-4*J 

PRINT 4ll,°EI,SPEI, XISUM,XISIG,KKK 

F ORHAK - l X-r*^^ *4 €4 2»4 r * STO£V=* , E12 .4 ~ r^-T I MEAN=* , E12.4 
+* I STD OE4=* ,£12.4»* ,9ITS=* ,E12.4) 


001602 


001604 

-004645- 

001637 

-641-644- 

001615 

.■004 64 .4- 

001623 


001666- 

001537 

00 16 4 5 - 

001645 

604645- 

001653 

0 0165 3- 

001663 


A 0 0 S= 3 . 


iMaafaiff>)7 


00 4G3 1=1, IDO 

AQO -A= A Q0 A- 4- -AQfl (-IV 

ADO S= AO OS +A DP (I)* AD 9(1) 


A0OA= AOOA/FIDO 

... A OQS^.SQgJ f.ADO S/E-XQQ-AQQ-A*-ADQA) 

PRINT 404, ILOCF, IDO, AOOA, ADOS 

F0RMA T (-4^,J4,-4SK 4N = * , I &,-* - , OROPWT-S^S46 , > , A VE = * , El 2 .4 , 

+* , ST 0 DEV=*,£12.4) 

_ 0440446 

PRINT 407, ILOCF 

FORMAT ( IX ,-*L-QCKlN=*- T I6-r*-rQRQPW TS- > 0*) 

CONTINUE 

— PRINT- 224 ,£• RR-M ~ 

FORMAT ( IX, *ERRM*,/, 10 (IX, 10E1 2.4,/) ) 

_f«^«:-343-,C-94N4,£ PENS — j 

FORMAT (IX, *-1,-1 P(6/A)TIMES*,F6.0,37X,*“1,+1 P (B/ A) TIMES*, F5 . 0) 


V-8 


00166? 


001710 


00 315 J=l,8 





315 F2 (J)=SBlNS2tI, J) 


316 FORMAT(iX,8F6.0,6X,8F6.0> 


001720 317 FOR|mATC1X,*+1,-1 PC 8/ A) TINES* ,F6. 0 , 37X,*+1 , U P(8/A)TIHES*,F6.0) 


001722 


DO 319 J=l» 8 



001727 


lIlbVAfcfl 


01745 


iVtifr/AI 


001761 


liiri.Tl 


01771 


iKniMil 


32000 


IfVfnhViltl 


002004 


319 F2 ( J) *S8I NS4 Cl , J) 


PRINT 304,11,12,111,122 




PKAV=PKAV/KDATA 


I it :4s% i FwHd ■ < li :« k TA :4 1 i u :<#< J VL2 U :4 • V k 1 


PRINT 305, PKA V,PKS0 




410 CONTINUE 


IJrJIJlflilfF 


104 F0RMATClHl,/7//> 



000006 
.t ,1 n n o <c 

o y u y 


DATA T , IRT/ 4194 304 • ,2 04 8/ 


000014 
■ 0000 30 

0 00 0 26 

non. -?a 
UUUUUl 


IYN=M0D<IYN,IRT) 
4 ?N^0 WI-B * tA-^ZN^I RT+- 


ROM^FLOAT <IYN*IP.T+IZN)/T 

3 T T I IQ M 

“ r. • r» 


000032 


END 





SUBROUTINE 91FILT (X ,Y »OHC ,T ,INIT) 

G INPUT INTO XC 4 ) , OUTPUT FROM Y4 4 ) , INIT =8 FOR NEW OMC OR T 

C FIRST OROER INTERPOLATION 

ww cmawMW 

000010 IFCINIT.NE.UX GO TO 1 

OO fl Oll INIT = 1 

000011 OHT=OMC*T 

000012 91= (( 0 rtT»8 . )*0MT/2 4 . »1 . X» 0 HTH . 

0000 20 92= (I (li.*0MT/24.*l.)*0MT-l.) +0HT-3.)/81 

■8899.” 8i=ltm , TQMT/2 4. wl .) T P HTm l, ) tQHT»3 . ) /8 1 

000035 B4=U(0MT-8.» *0MT/24.+i.)*0MT-i.)/9i 

000044 A l sQHT»0MT»0MT/<2 4. +Bl) 

000046 A2=li.*Ai 

C OU TPUT OMT , A1 , A 2y Bl »B 2 » 8 3,B4 

0 00050 1 Y(4) = A1*(X(4) +X(1) ) +A2* (X(3) + X(2) ) -B2*Y <31 -B3*Y ( 2) -B4*V (L> 

000125 nn 2 Tsl.,.3 

000126 X (I ) = X (I + i) 

4X410x34 2 rm = YCT»lT 

000144 RETURN 

000144 EM- — 


00 0 0 1 4 

000014 

000014 

000015 

000017 

000020 

000 0 20 

000022 

■000 * 33 

000025 

00 0 0 4 3 

000046 

0 0 0052 IX 

000064 

00006 6 

000074 

-0000- 75 10- 

000077 20 

00 0 -101 

000102 12 

0001 1 1 

000115 

0 00 116 

000123 13 

-00 0 1 30 

000132 

0 00 1 33 

000136 

0001 4 2 

000142 30 

-0 00 14 3 

000144 31 

- 08 0 4- 4 5 


SUBROUTINE 9LOBA(N,GAM, VAL,CCL, IL,S,9,IYN,IZN) 

D IMENSION TG4128) , RC<128) 

DATA G43/. 892979511/ 

IFm i.4 Q . 0 .» GO. I Q - 3* 

IF (CCL.EQ.1 •> GO TO 10 

EX=l . /3 . 

YN=0. 

S - * =-0-» 

00 11 1=1, N 

tf~FLQAT 4 1 ■! XT GAM 

TG(I)=GAH/G43*V**EX»EXPt-V) 

S 2 =S2*TG(I)*TG4I ) 

CALL NRANDM4IYN, IZN ,RC II) > 

YN=yN»TGCI)»P.C(I) 

S22=SQRT(S2> 

YN-YN/S22T UAL 

GO TO 20 

YN — YT 

NN=N-1 

- 0 0 12 J=1 , N N 

RC (N+l-J) =RC( N-J) 

CALL NRANQHdYN, IZN , R C| 1 1 ) 

YT=0 . 

4 X 0-13 JJ = l r N 

YT=YT+TG( JJ»*RC< JJ> 

YTtYT/S22»VAL 

FILY1=FL0AT (IL -1) 

S= 4 VT - YN)/FXLM 

B=-S*FLOATCIL> *YT 

GO TO 31 

S=0. 

-9 =0 - . 

RETURN 

-&WO 




SUBROUTINE BL0BP<N,GAM,VAL,CCL,IL,S,8,IYN,IZN) 


100014 


DATA G43/. 892979511/ 

VUUU1T 

000015 
n n n n i 7 


IFtCCL.EQ.l.l GO TO 10 

r V- 4 4 1 .. 

000 0 20 


cA-li / - — 

YN=C. 

oo— r\ — . - - 

nrv k)U cU 

000022 
rt a n a ^ 7 


> £ — v • — 

00 11 1 = 1 , N 

if » c ■ n > r f r a % » r A u 

V V w V u V 

000025 

n a n n /, ? 


TG(I)=GAM/G43*V**£X*£XP t-V> 

fl.tiiTCM* *TffTI .... . . 

-tf U U g^rO 

300046 

ri n n u s ? 

1 1 

3 C ~ tg T ±~T • J -4 • mil 

CALL NRANOM CIYN,IZN,RCd>) 

V»-V»UTf. f T\ *pr ( T> - - - - 

trv U U C 

000064 

1 X 

S22=SGRT { S2) 

VKI- VM /C>9*U Al 

000074 
fltinn 75 


GO TO Zii 

vm-YT - . 

U V WU f 3 

000077 

IV 

20 

NN= N- 1 

no 19 i -4 . uu 

000102 
.i ft n 4 4 4 

12 

— L iW lx Wt- 

RC<N+l-J)=RClN-J) 

( mo ninu it vm - tjm ormi ... 

000115 

A n ,'i 4 4 C 


YT = t. 

n n a ■? i i * a ki — 

tnnri rP 

00ul23 
n n 1 1 7 r, 

13 

1't 3 J J* 1 y n 

YT=YT*TG<JJ)*RC(JJ> 

VT*VT / C99*lfAI .. . ...... 

wwlvt 

000132 

nflm 1 ~z 


t t-H / J t i Vnl. 

FILM=FLOAT (IL-1) 

C - fwT Vkl\ ICTl M ..... - - - 

300136 

n .i n 4 /. o 


B=-S*FLOAT<IL) +YT 

no TO T1 - . . . .. 

UyU iHc 

000142 
"i n n i fi 7 

30 

S = 3 . 

n — n . 

000144 

31 

RETURN 

cun 

THU * TV L- 1 ~ "J 

J3000 7 


SUBROUTINE TRIG (CCLT, ANGR,S ,C) 

9IMFNSI0N SINEC91) - 

300007 


0R=. 0174532925 

TP f m T.f n. l 1 CO TO IQ - . 

a 3 0 012 

.innnii 


NT =91 

00 9 1=1 . MT 

■SvnriftEMPi 

9 

10 

SINE(I)=$IN(FLOAT<I- 1 >*OR) 
TP^XnT t AN^o /no *■ ,5 > 

3 00 0 32 

n ii fl a tc. 


IM= A9S ( ANGR/DR* . 5 ) 
T F (I M tGT. 2 7 0> G n T n 4 

■TlrTTiT^Hi 


IFdN.GT.18Q) GO TO 3 

IF(TN.GT.90) GO T 0 2 

000350 
414-04-51 


S=SINE(IH+1) 

3 0 J 053 
0 000 S ' 1 


GO TO 33 



000055 

u0Q05 7 


C= -SINE ( I M- 94 * 1 ) 
Gn tn j(] 

000057 
J 00361 

3 

S = -SINE(IM-18 0+1) 
C‘ , -STNE(-IM+2 704l) 

000063 

300053 

If 

GO TO 30 

5 = »<iTMf (-IM*T60 + 1 ) 

030065 

003067 

30 

C = SINE(IN-27C+1) 
IF(TO.LT.O) S=-S 

0 0 u 072 
4-04141X3 


RETURN V ’ 









SUBROUTINE A$ CP.TF ( X U, S) 

O 0 Q Q O^ - -X= ABS4-XU V 

0 3u 0 06 IF(X.GT.12*.> GO TO 1 

034442- IMX^-GT.64,*) . £0 TO- 2- 

03130155 TR(X.GT.32.) GO TO 3 

3 30-6 20--— — 1FIX.GT.16.) GO TO 4 

303023 I F ( X • GT.6.) GO TO 5 

0000 2-6 IF tX . GT .4 , ) GO TO 6 

000031 IFCX.GT.2.) GO TO 7 

6EU034 IF (X,GT*l-a--GJ3 TO & - 

300037 IFCX.GT..f)GO TO 9 

0 0 0942 -•• IF ( X* GT . .25) GO TO 10 

3030+5 IF (X . GT . . 125) GO TO 11 

444454- - s =-.-36 3*; 

300051 RETURN 

503454 i — G-i-11* 113’ . .. 

3 0 0 0 62 RET \V>H 

0 00055—2 6 = -. 0517, MX-128. ) *11. 313? 

300057 RETURN 

0*3057- 5 **,57-322* (X-5-+.)*d. 

3 00 063 R £ Tll° N 

466463 U - < v=-. 10455* (X=52. > *5.6564 

030067 RETURN 

0 00467 5--- Si. 14644* (X-1&.) *4, 

J 00073 RETURN 

0 0 0075— 6 S= , ?0?1* (X-5* ) +2.8264 

000077 RETURN 

6004-7-7 7- S^242 8* 4X-— ) *2 . 

000103 RETURN 

306103 -* S = .4142*(X-2. M-1.4142 

003107 RETURN 

003147 -9- • — -3*-, 6457* (X-l.) 4-1. 

0 00 112 RE TUP M 

0 SO 143 40 -.-826*+* (X-,4) +-.7 07 1 

* Qu 117 RETURN 

344117 —41 S-l»l7157*(X-.25) +.5 

000123 RETURN 

000123- - 


SU^F^UTINE OATASCIXU) 

004444- 44144491 ON IRC (6) 

003003 RATA T^C/w » v » C» 0 » u » 1/ 

000 0 0 3 4*--C- 

300003 IEfTRC (1) .NE. IRG (6) )IX=1 

3-64*-&7 tKV 1 -T*l,9 

OuOJll 1 TRC (7-1) =IRC(6-T) 

404B4X I4C <D = I* 

00002C IXUsl 

3000 2 1 IF(TX . E Q.464*U=-4 

000020 RETURN 

4 444 24 tW 


V-12 




COHT. ON PD 2 


PQ 21 


V-13 





















PH16S=PHtOEG«Pt/lBO. 
PH1SSB=PHIB«PI/10O. 
RMlN=S<JRT{ENOf1«NO>BRflTE ) 
AN IN2=AM IN/ 1 -4142 


At>C=SQRT( AMIN2 ) 

AMIN3=AMtN2/SQRTC 10 .•»( OYNR/ 10 . ) ) 
A=SQRTtEN0»NO»BRATE ) 


CALCULATE LOOP PARAMETERS C TRU1 . T AU2 .K . TAUB 1 .TRUB2 .KB ) 


2 

MN=0L2/tZETA+O.25/ZETfl) 

AK=P12«F0IFF/SINCPHISS) 

TAU1=AK/(UN«WN) 

TAU2=2.»CZETA-0.5/(WN«TAUl) J/AN 

KrAK/ANIN 

E00P=PI2«FRQ00P/K 

E00PR=P12«FQ0PRT-DT/K 

EOOPRXsEOOP+EDOPR 


EVnAXrPl2-SAMX/K 

EVniN=-PI2«SUtlN/K 

SVR=P12*SVRHZ/K 

UN8=BL2B/f ZETRB-0.2S/ZETRB 1 

AK=PI2«FD1FFB/SIN(PHISS9) 

T AUB 1=AK/ I WNB«WNB ) 

TRUB2 = 2 - •(ZETRB-0 -5/C WNB«TAUB 1 > 1/HNB 
RKBMIN=PI2*£0£LF/2- )/(K«l .4 125) 


K0=AK/AKBH IN 


C DEFINE CONSTANTS 


cu 

l.-DT/Tfllll 

C2r 

TAU2/TAU1 

C3 = 

C 1 »C2- 1 • 

C4 = 

OT-K 

cs= 

0T/TAU3 

C6=; 

1 --C5 

DTH 

I=Pl-OELF*OT 

| SOTHOrSlNtOTHl ) | 


C0THUC0SIDTH1) 
T)THX = PI2*FLO*OT 
SOThXsS INC OTHX ) 
COTHX=COS( OTHX 1 
C0l = l --BPER/TAUB1 
CB2=TAUB2/TAUB1 
C03=CB 1 ♦CB2- 1 . 
CSMS=0T/TAU4 


C6MS=l .-csns 
C5MM=DT/TBU5 
C6MH=l .-C5HH 
C50P=0T/TAU6 
C60P=1 --CSDP 
6IO=0.5-SQRTtNO/OT ) 


5 













INITIALIZE LOOPS (tH$.EO*XI>XQ,QV.QERRO.ERR) 


6LSUfl-0 ♦ 

SLSIO=0 . 

XISUP1 = 0 . 

X IS 1 0-0 - 

NB I T = 8PER/ 0 T ♦ . 5 

NBIT2=N0lT/2 

TSYNC=0PER 

CAMPC I ) = ( 0 . .0 . ) 

Tr~ z . 

C AHF(2 3--I 0 . .0 ■ ) 

c ahp 1 3 1 = c a - .o . ) 
CFILTC i ) = (0 - -a ■ I 
CFILH2) = (0..0- ) 
CFILTC3) = lO..O - J 
I N 1 T = 0 

OflCrPI-BiF 

IOflTfls-1 


PRINT 201 -EN0D8.0ELF .B IF .012 -NS LICE .AS OFT 


r 0RH A TUX . *ErOC3= ■ .F5 . I 
• * BL 2= « . F S .0 »• .N5LICE 


c LF = * . F 5 ■ l . ■ .31F=«.F6*0« 
, * . ASOF T z ■ . F S .4 ) 


PRINT 202 -SVRHZ .S WnX .SUriN .FLO .FR JOGP .HJOPRT 


FORnATt IX . ■ S VRHZ = • .E12 .4 .■ . SWnX = ■ . E 1 2 . 4 . ■ .SWrtNz-.Et2.-1 
■ .F|_0 = « .E12 .4 .• ,FRQOOP=. .£12 .4 .« .FQOPR T = ■ .£ 1 2 • 4 ) 


PRINT 203 .TL .QVNR .SORnP .SOPHfl .11. 12. III. 122 


FORnflTl IX .«TL=».F6 -2 . • . D T nR = » . F6 - l . 

■ ,SOflnP = « ,F6 .3 .» .SOPHfls. .F6 .3 -■ . I I . 12 . 1 U . 1 22« - 4 16 ) 


PRINT 204 .UN .K .TAUl . TflU2 .UNB .KB .TAU31 . TRU32 


FORnflTl IX .»UN = . .F6 • 1 .• .« = ■ . F 6 . 1 . » . TflU Is. .E 1 2 ■ 4 , • . TRU2 = ■ -E 1 2 - 4 . 

- .«N0 = -.F6- 1 .• .K0I..F6-1." . TflU8 l = ■ .E 12 ■ 4 . ■ .TflU02=- ,E12 .4 ) 


IERRSP=0 


CONT. ON PO 














































I 



EP-7 OF g3 


V-19 
























































V - 2 3 








ro 














XIO = X1/FLORTCN0IT I 
vnH=(CGnH»VMHfC5rtH«unH J-flS 
uuhsO. ■ 



























V-27 























B5 




CONT. ON PC 18 



V- 29 








corn • on po 


19 




































































I 




V-36 





SUBROUTINE Rf)NO(1( IYN.IZN.ROfl) 
OftTfl Ifl . IB/129 . 1/ 

ORTA T.IRT/4194304 . .2048/ 
IYN=IflnIYN*[ lB*Ifl«IZN)/IRT 
IYN=H00t IYN.1RT) 

IZNsflOOC [ IB* Ifl« IZN ) . IRT ) 
ROHsFLOAT ( I YN« IRT* 1ZN I/T 


ITeturn I 



SUBROUTINE eiFllTtX.Y.OnC.T.INIT) 



PQ_I E-LMRL 


V-38 




























































V-42 











V-43 
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V-51 









APPENDIX VI 
ACQUISITION 
COMPUTER SOFTWARE 


VI-1 



3 00003 
000003 
300003 
000003 
3 0 0 0 0 3 
300003 
000003 
300003 
0 00003 
C 0 0 0 0 3 
300003 
000003 
000003 
000303 


00 0003 
DQ0006 
000010 
0 00013 
000016 
0 000 20 

0 C 0 0 2 n 
300023 
Cfl 00 24 
00002* 
Q 0 0 0 26 
003010 
000031 
0 013 0 32 
J0003T 
000034 
3 0 0 0 36 
000037 
000040 

0 0 PO 4 ? 
000060 
000066 
J 00 060 
300062 
330064 

ouoan 

OCOO-^3 

000100 

003102 

000113 

003117 

300123 

000126 

000132 

0001 34 

VI-2 


opOG^AM MfiRKfTNPUT, CUTPL T ,T A 0 E6*IMPUT ,T APE 6= OUTPUT) 

REAL NS, Nr, NO ,K,KP 

tOMPl^X CAMPT4) ,CFILT (4) ,fOUT ,CrOUT 

DATA *3AMP,FRPHA,SOAMP, EOPHA, N0/2.,.2».23,.47,l./ 

DATA GAM43,N43/.09,64/ 

DATA ENfJ ''E, 6P ATE,EN CM IN/ 3 . , 80 . » 6 . / 

DATA DElF,NSPB,NAC0T»6IP/62. ,40, 100*10 CO. / 

DATA «L2, PHT n FG, FDI FF/176 . » 10 . » 34 2. 433/ 

DATA T1,T2, 111, 122/1,1,1,1/ 

DATA FRQDCP ,FOD°RT,TL/i3.E3,10.8, -6 .5/ 

OATA TAD3/ .07 967/ 

DA t A PLD/10./ 

DATA 6VRH7,SWMX,SWMN/2. 2E3, 28.16E3,6.6FE3/ 

DATA DYMR/3 • / 

NAHELI?T/IN1/MAC0T, NSP8»EN0DR,i5ELF, OIF, EMCHIN 
♦,6L2,°HT0EG,EniFF 

♦ ,F^AMp t FiPH A, SDAMF,SDPHA, PAH43 ,N43 
+, FLO, FRQDCP,FnO°PT, 11,12, Til, 122 
+ , P W°H 7 , SWMX ,S WMN , TL ,DYNP 
PT=L. *ATAN(1.3 
P I 2=2 • #P I 
t ft a=pqp,t ro .= ) 

rftU TRIG (CrLT,°I,SSCIN,C6riN) 

CCLTsi. 
r ONTI NUP 
REAP INPU t S 
RFAP(6,T K U) 

NSPP* 200 
M AC C T -3 0 

PWMX=260Q. 

connoP = 100!3 . 

T 1 = 1 
12=1 
111=1 
T 2 2= 1 
T1S=T1 

I 2 S= I 2 
T 1 1 C =T11 
T 22*= T 22 

CONVERT INPUTS CEN0 ,BPEP, PHTSS ,PHTS66) 

F N 0 = 1 C.** (F M0 DP/10. ) 

ENC«=10 .** (EM0MTN/1C.) 

D r FP= 1 * / DP A T E 
D T = FPE 0 /MSP n 
ohIS6=phtdEG*°I/18C . 

AM IMPORT (ENP M*M0*PRATE> 

AMTN2=AM t N/1,h142 
A = 5 QRT (E M0* M0 * c * p ATE) 

AG C=6 CRT C AHTM2> 

AMTM3=AMTN2/GDRT(10 { ny mo / 1 jj . > ) 

C ALC'IL A T E LOOP PARA METERS ( TAU1 , TAU2 , K , TA UB 1 , T AU° 2 , K °) 
WN=«L 2/( 7ET A* 0 • 26 /ZET A) 

AK = FI?* C DT PF /SIN(PHTS6) 

TAU1=AX/(WNTWN> 

TAU2=2.»C ZETA-0.5/tWN*T AU1) ) / WM 

K- AK / A Ml M 

E p 0P = PI2*FR0n OP/K 



300137 

£D0PP=PI2*FQ9RRT*0T/K 

990142 

€VWAX=PI2*SHWX/K 

000144 

EVMIN--PT2* SWMN/K 

300146 

SVR=ei2*SVRHZ/K 

JLIU11S.0._. 

DES=9T*SVJ? 

c CEFTNE CONSTANTS 

000152 

Cl«l» -CT/TAIJl 

090155 

C2=TAU2/T AUl 

000156 

C3*ri»C2-l, 

909161 

C4=OT*K 

900163 

C5=0T/TAU3 

300164 

C6=l.-C5 

000165 

OTHi=°T*OELF*DT 

900170 

SOT HO = SIN (OTH 1) 

000172 

COT Hl=COS (OTHi ) 

300174 

OTHX=PI2*FL 0*DT 

00017T 

50THX=STM(0THX) 

090201 

COTHX=COS (OTHX ) 

300203 

NAC 0=0 

000204 

ACO-O * 

000205 

A VEEF = 3 , 

U 00 206 

STOFS=fl. ' 



0Q0207 

aoozic 

0 0 3 211 
000 212 
300213 
000221 
.0.00234 

000234 

000254 

000254 

JJ00330 

000300 

000314 


.000314 

000314 

00032? 

900322 

000323 

000347 

J1M347. 

300347 
000351 
000354 
000356 
300361 
HO 036 3 
300366 
000 370 
0 00 3 7 3 
300373 
003375 


A^OU=0. 

SKOU=0. ’ 

AMISS=0. 

AFAL = 0 • 

SIG=0.5*SGRTTN0/DT) 

°RINT 202 ? EN0O'R,DELF, EIF,RL2 

202 FORMAT (IX, *EN0OR=*,F6,l»» , DELF=*,P6. 1, 

♦* ,RIP=*,F6.0,* ,9L2-*,F6.0> 

PRINT 204 ,SVRH7,S HMX, SWHN ,FLO , FRQOOP, FQDPRT 
204 FORMAT fix ,*5VRH7-*,E12«4,* ,SHMX=*, E12.4t * , SWMN=*, E12. 4, 
,FL0=*,E12.4»* ,FRQ00P=*,E12.4,* ,FQDPRTs*,E12.4) 

PRINT 206 ,TL» DYNR , SOA HR , SOP HA , I IS , I 2S, IliS, I22S 

206 FORHAT(1X,ttl=*,F 6,2,* ,HY.NR=*,F6 .1,* t S&AMP=*,F6 .3, 

+* »SOPHA=»,F6.3»* , Ii, 12, 111, 122*, 416) 

PRINT 2(J7 , WN, K,TAUi jT AU2 

207 FORMAT(iX,*WM*»,F6.i f * ,K=*,F6.1, 

♦* ,TAUi=*,E12.4,* »TAU2-*,E12*4) 

- INITTALT7E LOOPS CTH2 ,E0 , XI ,XQ, QV, QERR0 ,ERR> 

4QD CONTINUE 

IF(NACG.EQ.l) PRINT399 

799 F9RM&T(lX,*NACQ,AMlSS,AFAL,KOATA,ES,EV,EVi*> 

IF(NACQ.EO.O) G0T0396 

PRINT 797 f NAC C, AMISS, AFAL *ACQ ,KDATA ,ES, EV ,EV1 
T97 FORMAT (IX ,T6,3E12»4»I6»3£i2»4) 

396 CONTINUE . 

NACQ= NACO+1 

IF (NACQ.GT.NACQT) GOTO 40 2 

CANPf 13sf 0 .,0 •> 

CANP(2>*(0.,0.) 

CAMR{3) = (0.,C .) J 

.CFILT (13= (0.,JJ.) 

CFILT(21-(0»,0») 

CFILT (3)* (0.,g.) 

INIT=0 
OMC=RI*0IF 
TO ATA=-1 


VI-3 



3903?6 


ll-Q . 

J G 0 3^7 


F V 9 = 0 . 

300430 


"H? = 0 . 

J004G1 


PTH2=1. 

103433 


5TH2=0. 

300433 


°THy= l. 

3 0 3 4 04 


<5THy= o • 

0 0Q4Q5 


5 0 — P • 

QC94S5 


EV=C. 

0004*7 


£onp^y=Enci i »5nnpR 

10u411 


EVl=E v+eoofpx 

300413 


51 = 1. 

000414 


G?*i. 

300415 


V= • 1 

3 0 341'* 


E F = C . 

3*041? 


5TH1=C , 

000420 


F ?H1=1. 

■3 0 34 2? 


Jr C 


r 

TIIPPULEN^F INITIALIZATION 

33*42? 


ACOP=SCR T (2. **.75-1 .) 

320432 


AAMP=Pl2*F n AM Pf AC OR 

3004 35 


Ar^)^ = DX2*copM A^ACOR 

3 0 04 3"’ 


TAMF=IN?( ,5<-G4M43/c AAHP*Q t ) ) 

3 30444 


TPH4=TNT( ,5*GAH43/( APFA*OT) ) 

00*44? 


r,AM43A = n T *FtnATtIAMP) *AAMP 

300452 


G ft M4 3 p=DT *FL 0 A? ( TPH A) * A PH A 

0 Q *>+5 e 


531=0 . 

30045"’ 


‘■'ALL PLC3AfN4 3,GAM4 3A,SrAMP,CCL,IAMP f 5lOA,8AHP,Ill,I22) 

00045" 


C Alt 2tn3F(N4 3,GAM4 3P,SCPHA,Ca*IPHa,SL0P t BPHA, 111,122) 

o o o F :• o 


cn»i . 

9905*1 


r AMrr = g 

300532 


I°HAO=9 


r 

5E1F5? DATA f IOATA) 

3005"? 


K p A ? ft = 0 

9 3 " 5 9 r 

4 

PQNTTK'UE 

303535 


IAMPP = I 

3 33 r 37 


T°HftC-I r H AC ♦I 

] 305 l r 


TF(IAHPC.GT.TAMp) G 0 TO 10 

3CQ513 

13 

IFCTPHAC.G? .TPHA) GO TO 11 

30051"’ 


GO TO 12 

00351? 

10 

T AMp^l 

3 005 25 


* All 91 09 A <N47, G AH43A,SOARP,CCl»I AMP,SLOA,8AMP,Ili, 122) 

030531 


GO 13 

300532 

11 

I °H AC=1 

330533 


C A 1 1 81 085 (H4 3, GAH4 3»,SPT>HA ,CCL » I P H A, SLOP , 8°H A, 111 , 122) 

300544 

12 

rON?Tft'l)E 

300544 


VA = S10A*F10 A? (IAHPC) ft-PAMP 

190547 


V C = S1 C°»F10AT < I PH AC ) * Q P HA 

0UC5 C ? 


■'flu ?RTG (O^t T,YP,SSCIW»C5CIN) 

300556 


TFfJ.lT.NSP8) GO TO 1 

3 C 0561 


J = Q 

0 1P561 


CALI OA T A5( IOATA) 

30056? 


SnTRi = SDTH0*TT)ATA 

0 P 0 56 5 


KO A?A = KO A T A *1 


<■> 

O £F TNF TNPUT5 ( T H1,NS,NC) 

3 0 C5&" 

1 

’'ON?? ft'UE 

33356? 


T FHF=S T HX*COTHXfCTHy*SOTHy 

3 CO 57 ? 
VI-4 


0THy=CTwy*C8 T Hy-STH x* 5ri th x 



903575 

900576 

030602 

000604 

0-11 0634 

00061b 

000625 
0 0 06 27 
900631 
00063b 

3JU3S_3_e_ _ 

000640 

000657 

000662 

000670 

L0 0672 _ 

000676 

000676 

000^03 

000706 

000710 

000711 

000716 

100723 

000730 

000^32 

000736 


030740 

000^41 

000744 

000750 

000753 

000764 

a0!P56 

000757 

000764 

000765 

000765 
90077Q 
000772 
000773 
030776 
001000 
001092 
001093 
"91005 
001006 
Q 0 10 13 
001315 
001022 
301024 
001032 
001041 
001041 


STHX-TE MP 

TEMP=STHi*CDTHi+CTHi*SnTHl 

CTHi=CTHi*COTHl-STHl*SDTHl 

STHlsTEHP 

AS_T Hi=l . 4142* A* C1.+YA1 *J STHJ*.C SCIN*£TH1*SSCIN1__. 

ACT HI =i, 4142* ft* (i«*YA1* <CTH1*CSCIN-STH1*SSCIN> 

CALL NRANOM (I1>I2»X1 
NrstX*5TG 

CALI NRANOH (Tl t I2*X1 
NS-X*SIG 

J-J*l . 

r COOP EQUATIONS 

CAHP (4) =CPPLX( ft STH1#NC, -ACTH1-NS1 *CNPIX ICTH2, -STH2) *G1 
CALL 91FILT trftMP,CFTLT, CMC,DT ,TNIT> 

CC0UT = G2*CFILT(4) 

OOUT-=CCOIJT 

V-_C6*V+C5*U 

VD = V 

IF ( ARS( VD1 .LT.AMIN3)V0= AMIN3 
IF ( VO* N£ • 0. 1C ALL ASQRTF (VC, VSQRT) 

Gl= AGC/V5QPT 
G2=Oi 

U=ATMAG(CCUT1 *C T HX-REAL (CCUT) *STHX 

E=PEAL ( CO U T 1 *CT H X+ A IM AG (CCUT) *STHX+ES 
E V=C1*EV+C2*E-C3*E0 
EVl=EV0-EDOPPX 
IF (FV »GT. EVMft XIGOTO500 
IFCV.LT.TH GDTO401 
C _ IF(EV .GT.EVMAX)OES*-OES 
r IF(EV.LT.EVHTN1T5ES=-fl€S 
EOOPRX=EOOPRV*ECOPR 
TH2=TH2+C4*EV1 
TH2=ANOO (TH2, PI 2) 

CALL TPTG(CCLT,TH2,STH2,CTH2> 

EVQ =EV 
EC=E 

EV1P=EV1 

IF (V.GT.TL>E5=E54- DES 
GOT 04 

401 OONTIKUE 

C HAVRE AN ACQUISITION 

TFTES.LT.EtJOPRX) GOT 0501 

ACQ= ACQ + 1 

AVEE*= A VFES*ES 

STOES*STOESbE S*ES 

AKOU- AKOU*FLO AT CKOATA) 

SKOU=$KOU*FLOAT(KOATA)*FLOAT(KOATA1 

GOT 04 00 

402 PAOO= ACQ/FLO»T(NACQT) 

IFTACQ.EO.O.) GOTO504 

STBEV=SQRT(»ACQ*Ci.-PACQi /FLOAT (NAC QT) ) 

AVEES=AVEE5/AC0 

STOES=EQRT(STOES/AC G-A VEE5*AYEE51 
AKOU= AKOU/ACQ 

SKOU^SQRT (SKOU/ACQ-AKOU*A«OU> 

PRINT 4D3,PACQ, STDE V 

403 FOPHAT(lX,*PACQ=*,E12.4,* , STDEV=* , E12. 4) 

PRINT 404 , A VEES, STDES 



QOiGSi 

QQim 
Q01061 
901061 
101061 
°9l364 
J 0 1 071 
0010^1 
oaio ^ 
001101 
001101 
00 lil c 
031115 
001121 
301121 
001122 

001122 
001124 
0 n 1 1 2 e > 

00112 r 
30112^ 
031130 
001134 
00 1174 

00 113 r 


0009J6 

300037 

□ □ a c i r 

0 0 0 G 1 1 
000017 
030017 


330006 
0 30036 
0 Q0006 

300314 

03302° 
•3 C 3 0 26 
3000 71 
J Q 0 0 '*2 


404 FORMfiTfix t MVEES-*, E12.4,* » STP FS=* , El 2 .4 ) 

PPIMT 40 c , AKTUJ, SKOU 

405 FORMAT(iY,*AVE KOAT A=*,E12.4,* »STQ OEV=* ,E 1 2 . 4) 
506 PONTTMJF 

PMTS^-AMICS/FLOAT (M ACQT) 

°P. T NT 50 2,PN T SS 

502 ropMAT ay ,*PP MISLED ACG=*,E12.4) 
oFftt-AFAL/FLOAT IN ACQT) 

P^INT 503 , c FAL 

^03 FORMAT ( IX ) fDO FALSE ACQ=*,E12.4> 

ORTt'7 304,11,12,111,122 
304 format (IX, * 11, 12, Til, 1 2 2*, 4 16) 

307 OD IN7 104 
U4 FORMAT (1H13 
GO Tn 3 
500 CONTIKUF 
C MI^S fP ACO 

AMiso = fl?iTCC<fi . 

G0TO4 00 
5 01 r OMTT NUE 
r r ALEE ACQ 

AFAL= PFAL +1 . 

G r >TP4 G 0 

504 D PINT 5 00 

505 FORMAT (IX ,*NO ACQUISITIONS* ) 

GOTO c 06 

EMO 


5U6 p 0U T TNE NP ANOMfl YN ,1 2N,ROMN) 
fnMN =-6 . 

Q O l N = 1 , 12 

P ALL PANQMdVN,TZN,ROM) 

1 RnMM = POMN+ROM 
RETURN 
F M 0 


ntPPOUTTNE RA MOM ( IY N, I 7K , ^OM) 
OATA IA, 10/120,1/ 

O A T A T,T°t/41 94304. ,2 048/ 
IYN=T A*IYN*(T9*TA*I7N)/IR t 
IY N =« 0 O ( T YM , tpt ) 

T7N=mPQ< (T9mTA*T7N) , IRT ) 
RDM=FLOAT <IYN*IRT*T7N)/T 

return 

FNR 


VI-6 



SUBROUTINE 91FILT(X,Y,0MC,T,INIT) 
r INPUT INTO XH4), OUTPUT FROM YI4), INIT*0 FOR NEW OMC OR T 

t FIRST CPH6R I NT ERPOt ATI CN 

300010 COMPLEX X(4),Y(4) . . - - 

030010 IFCINIT.NE.O) GO TO 1 

300011 INTTri 

000011 OMT*ONC*T 

300012 31= 1 1 0MT» 8* )* OHT /24 .*1.)*0MT*1. 

00 00 20 B2=T(til.*OMT/24.*l.)*OMT-l.)*CMT-3.) /B1 

Qfl(ia27 _ 33=C('CU.*0HT^24*jr±j).*0NI-l.l-?0MI»:l,)Zai._- — 

0 0 00 3 C B4= (((OMT-8.) *OMT/24. +1 . >*OMT-l .) /B1 

300044 A1 = 0M7*0NT*0HT/ C24.*B1) 

330046 A 2= 11 * * A1 

C OUTPUT 0MT,Al,A2,Bi,32 f 33,B4 

0 00050 1 Y t4)=Al* ( X(4) *X( 1) ) 4A 2MX C 3) *X ( 2) ) -B2*Y (3) -B3*Y 1 2) - 04* Y < 1 ) 

00 012* _ DO 2 .1=1,3 ... _ - - 

030126 X (I) = X<I»11 

000134 2 YfI)=YCT*i) 

000144 RETURN 

000144 ENO 


000014 


SUBROUTINE 9L0BAfN,GAM,VAL,CCL,TL,S,BtIYN,IZN) 
DIMENSION T Gf 128) *RCC128) 

0 C 00 14 


DATA G43/ .842979511/ 

10 0 014 . . 


IFCVAL.EQ.O .1..G.Q TO 30_ . .. 

33001* 


TFtTPL.EQ. 1.) GO TO 10 

003017 


EX=l./3. 

0000 20 


YN=0. 

000020 


S 2= 0 * 

000022 


DO 11 1=1, N 

300023 


V=FlOATfI-11*GAM .. 

000 0 25 


TG (T) =G AM/G43*Y**EX*EXP C-V) 

30004? 


S2 = S2*TGC I)*TG(T) 

000046 


CALL NRANOM (I YN , IZN ,RC (111 

000052 

11 

YN=YN+TG(I)*RC<I) 

D 00064 


S ?2=S QRT ( ?2) 

000066 


VN=YM/S22*VAL 

000074 


GO T n 20 

0000^* 

10 

YN=YT 

0 0 0 Q T 7 

20 

NN=M-1 

300101 


DO 12 J=1,NN 

00010? 

12 

RC(N*i-J) =RC(N-J) 

_0 00111 


CALL KR_ANOMTYN,IZN,RCCl) ) .. . 

000 115 


Y T - 0 • 

000116 


DO 13 J J= 1 , N 

00012? 

13 

YT=YT4TGCJJ)*RC(JJ) 

000130 


YT = YT/S2?*VAL 

00013? 


FTLM=FLOAT(IL-l) 

POP 133 


S= (YT -YN> /FIL M . .. _ . . _. 

000136 


0=-S* FLOATC II ) *YT 

000142 


GO TO 31 

000142 

30 

s=o. 

000143 


9=0. 

300144 

31 

RETURN 

000145 


END 


VI-7 



'J U 'J 1 4 


‘JintHMUM T G C 1 2 * ) ,RC\ 12 rr 

300014 


DATA G43/. 992979511/ 

'3 030 14 


IFrVL.EO.O.) GO TO 30 

31011* 


TF (C*! .EO, 1 .) GO TO 10 

30 3017 


EX = 1./3. 

0G0Q2C 


VM=G . 

000020 


S2=C. 

000022 


OO 11 T = l t M 

000127 


V=FLnfi T (I-1)*GAM 

3*00 2* 


TG (I) =GAM/FI*7*V**EX*EXP <-V) 

000147 


S2 = S2fTG<T) -fTG(I) 

J3u046 


CALL NRAMr^(TYN,I2N t RCtT) ) 

0 0 3 0 c 2 

11 

YNxYM+TG ( (I) 

] 03 364 


f 22=FRRT tS21 

i 0 3**6 


Y M=v N/F22* V AL 

130074 


r,o to jo 

1 0 00’’ F 

10 

YM=v» 

3 010” 

?0 

MMrt-i-i 

0 0 3 1 C 1 


OO 12 J=1,NM 

303132 

12 

c *'MMfl-J)=or(N- j) 

000111 


C ALL NR AN DM 11 YN » IZN jRF ( l ) ) 

’Gill* 


Y’ = 0 . 

000116 


^ 13 JJ=1,*' 

3 1 u 1 2 7 

13 

Y T =YT yTGC JJ)*RC ( JJ) 

013170 


YT=YT/F2?* VAL 

033132 


FTLM^FLOAT(TL-n 

000177 


S- (V T -YM) /FIL * 

3C*id* 


0= - S*FLGA t (IL) fYT 

3011*2 


GO 7^ 

a o c i u 2 

“»u 

f = 0. 

33*143 


«s0. 

*311*4 

31 

3 C- TIP N 

1 1 ' 1 1 45 


cun 


SURROL'TINF totgcCCLT, ANGR t S f C) 


03030’ 


DIMENSION FT^'E (91) 

30331’ 


D R= « 1 17 4* 32935 

?03n n 


T F (C r LT . r C • 1. ) GO to 10 

3 3 00 1’ 


MT= Qi 

3 13ul3 


DO 9 T-1,wt 

100014 

q 

F INF( T3-9IN(FL0AT( I-i ) *CR) 

10102’ 

1C 

I°~IN T C ANGR/^R*. 5) 

IIOO'I 32 


TM=ARF J ANGR/9R«-.5> 

00033* 


TF(I M .G T . 2701 GO TO 4 

000042 


t F(TM.GT. 160) GO T 0 3 

Jd 3045 


I*(TM.GT.90) GO to 2 

0001*1 


S = S I N E II M +11 

300051 


C = ftne ( “I **f 9G *-1) 

301053 


GO T 0 30 

003353 

? 

^-FTMEt-I^lF flf 1) 

31035* 


r = -STNE(TM-90fl) 

00005’ 


GO Tn 30 

301057 

T 

<?=-ctne(T m -1»0*1) 

0 3 9J61 


, '=-FjN E (_ I Mf270fl) 

J 3 J 163 


GO TO 30 

100053 

IjL 

S=-*IKE(-T*f360fl) 

000365 


f = ftiv)E (TM-2’Q*1) 

30306’ 

30 

TF(IR.I T „0> f = .c 

0 C 0 3 7 2 


OpTMON 

3 0 J 1’7 


CMP 
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€ 

SUBROUTINE ASGRTF (XU,S) 
MAX IN?ER P E?R=1»46 <? l.LE 

000095 


X=ABMXU) 

.0.0 0 CIO 6 


IFtX.GT.i28.) GO TO 1 

000012 


IFCX.GT.64.) GO TO 2 

000015 


IF(X»GT.32.1 GO TO 3 

OTJT30 20 


IFIX.GT.16.) GO TO 4 

000023 


IFIX.CT.8.) GO TO 5 

000026 


IFtX.GT.4.) GO TO 6 

_9_Q_Q_d 3 1 . 


IF(X,GT.2.1 GO TO ? 

300934 


IFtX.GT.i.) GO TO 8 

000037 


IF(V.GT..5) GO TO 9 

000042 


IFtX.GT..25)GO TO 10 

9OO04 K 


IF (X.GT..1251GO TO 11 

aooosn 


S= .3535 

990051 


RFTURN 

OOOQ51 

” 1 

S = ll. 311? 

000052 


RETURN 

000053 

2 

S=.05 1??*CX-128.) *11.313? 

000057 


RETURN 

00005T 

3 

S= • IJ7 3 22* (X— 64 • 1*8. 

000063 


RETU°N 

000063 

4 

S=. 10755* TX-32.) *5.6568 

000067 


RETURN 

000067 

5 

R=,ib644* (X— 16 • ) +4 • 

000073 


RETURN 

00907"* 

6 

S=. 237 1* t X— 8. ) *2.8284 

000077 


RETURN 

00807? 

7 

B= • 29 28* ( X-4. ) *2 • 

000103 


RETURN 

000103 

8 

S=.4142*tX-2.) *1,4142 

000197 


RETURN 

otaio? 

9 

S=, P 86 T *(X- 1,1*1, 

000112 


RETURN 

000117 

19 

S=.8284*(X-.5)*.707i 

90011? 


return 

000117 

11 

Ssl.lTlE’TfX-.aS) +.5 

000123 


return 

900123 


END 


*00003 

SUBROUTINE OSTAStlXU) 
Ol ME NS I ON TRCC61 

909003 

OATA IRC/C, G, 0,9,0, 1/ 

DO 0.0 0 3 

..... ix=o 

000003 

I F<IRC (1) .NE.IROC6) >IX=1 

000007 

OO 1 1=1,5 

000011 

1 IRC 17 -11 = IRE (6-1) 

00001? 

IRCT1 1 = IX 

000020 

IXU=1 

.aiCLUii— 

_ IF C T X . £0 . 0) tXU=-l . . 

000023 

RETURN 

000924 

ENO 



PROGRAM MARK! INPUT .OUTPUT .TAPES= INPUT . TAPE6=0UTPUT 
i REAL NS.NC.NO.K.KB 

COMPLEX CAMP 14 ) .CF ILT( 4 ) .COUT«CCOUT - . . 

DATA FBAHP .FBPHA .SDAMP .SOP HA .NO/ 2 . • *2 • -23 . -4T . 1 •/ 
ORTA GAM43.N43/ -09.64/ 

; DATA ENODB.BRATE.ENOMIN/8-.08*«9;/ , 

DATA OELF .NSPB .NACOT ,B IF/62 . .40.100.1000./ 

1 nflTfl BL2.PH10EG.F0IFF/176..1Q..342.433/ 

£ 

I DATA II. 12. III. 122/1. 1.1.1/ 

DATA FRODOP .FQDPRT .TL/ 13 .E3 . 10 .8 .-6 .5/ 

DATA TAU3/ .03957/ 
noTo pi n/ tn./ 

DATA SVRHZ.SWMX.SHMN/2.2E3.28-1SE3 .6 -8SE3/ 

DATA DYNR/3 . / 



I NAMELIST / INI /NACOT .NSPB .ENOOB .DELF .B IF .ENOMIN 
.BL2.PHI0EG.F01FF 

.FBAMP .FBPHA.SDBHP.S0PHA.GAM4 3.N43 
.FLO.FRQOOP.FQOPRT.il .12.111 .122 
.SVRHZ .SNMX.SNMN.TL.OTNR 


Pl=4 -»ATANl 1 . ) 
P12=2.»PI 
2ETA=SQRT (0-5 ) 


1 CALL TRI GICCLT.Pl . SSC I N . CSC IN ) j 
I 



CONT. ON PO 2 


PP 1 OF IZ 


VI-10 









CONT. ON PO 3 


PP 2 fl£ VL 


VI-11 






















VI-13 















CONT - OH PO 6 


ro S DF 


1Z 


VI-14 













5 




















CflrtPt 4 )=CHPLXt RSTH1-NC .-flCTHl-N6)" C rtPLXtCTh2 .-ST M2 )»& l_ 

^ — > 

fCRLL BlFILTlCftnP.CFtLT.OHC.OT.INtn j 


CC0UT=O2«CFtLTt4 ) 
COUTsCCOUT 
V=C6»V-»C5nU 
VO=V 


lF(fl6S( VO ] -LT .AMIN3 ) 


V0=flHIN3 


IFtVD.NE-O.) 


CALL ASQRTFC VO .VSQRT 




U=fllMHClCOUT )«CTHX-REflLtCOUT )»STHX 
E=REflLtCOUT )»CTHX»fl IMflOt COUT )»STHX+ES 
EVsCl»EV*C2»E-C3«E0 
EVUEVO-EOOPRX 


















rooTosoT :x® 


flCOrflCOtl 




flVEES=flVEES-ES 
STOE6=STOES+ES«ES 
RKOU=flK Oil ♦FLOAT t KDATA ) 
8KOU=SKOU«FLOAT ( KDATfl 1-FLOAT C KDATA ) 


2 


| IT0T0400 


flZ. 



CONT- ON PO 11 


Efl—ULflE 12. 


VI-19 














































VI-22 






SUBROUTINE RANOH( IYN.1ZN.RDH) 
OATH IA. 18/129.1/ 

DATA T.IRT/4194304..204B/ 

IYN = 1A« lYN-»( 10-IA-1ZN)/IRT 
IVN = nOO( 1YN.1RT ) 

1ZN=H00(( 1S+1A*IZN) . 1RT ) 
ROH=FLOAT ( IVN»1RT* 1ZN )/T 


Return | 




SUBROUTINE BlFlLTO .Y .OMC .T .1NIT ) 


INPUT INTO XI 4 ) . OUT'pUT FROM YM1. 1NIT=0 f 0R NEW OMC OR T 
FIRST ORDER INTERPOLATION^ 

1 COMPLEX XT 4) .Y(4 )1 


IF(INIT-NE.O) 

dr?EE> 

r 

i IN t T= l 1 


OMT=OMC«T , 

B1=((OMT»8.)«OMT/24.*1.)»OMT«1. 

82 = ( t ( tl-"0MT/24.+t . )«0MT- 1 • l-OMT-l . )/81 
B3=ri(ll.»0nT/24 - -1 .I-OMT-l . 1-0MTO- l/Bl 
84= J U0MT-8- J-OMT/24 • 1-0MT-1 - )/Bl 
flUOMT-OMT-OMT/(24 .*0U 

R2= 1 1 ■■fll 

y , 

ft OUTPUT OhT.Rl.fl2.Bl.8g.83.BM 


Y14 1sfil»lX14 1*Xm)«ft2"tXC3)Oit2) )-82»Yt3)-B3«Y(2>-B4«Yt [T 


: oo z 

^3 

1=1.3 > 
l 

1 xt n=xt mi 1 

2 < 

7 



RETURN 


i rera i 
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VI-35 






































APPENDIX VII 


MISCELLANEOUS SUBROUTINES 

This appendix presents three subroutines which were found useful in trouble- 
shooting the receiver during the software development, but were not used in the 
production runs. 

The first is OSCILL which is an oscilloscope and plots three inputs. These 
are entered into the routine by lists. The routine automatically scales the 
output. 

The second is SURV which is a voltmeter. It gives the mean and standard 
deviation of the input, as well as the mean and standard deviation of the absolute 
value of the input. 

The third is CHMEM. This was used early in the analysis when there was one 
conditional probability matrix which was the position by position summation of 
all four of the matrices of this report, i.e., just P{B/A), It generates infor- 
mation to determine the skew of the matrix. 
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SUBROUTINE OSCILL(LLL,X,Y,Z) 


C 

XYZ piotteo vs ut 

0 0 0007 — 


DIMENSION OUT (66) * X (90 1 ,Y(90> ,2(901 

000007 


OATA BLANK, STAR, XX, YY ,Z2,AZ/ 



♦ IH ',iP*,lFX;iHY,iH2,lHO/ 

000007 


L=LLL -1 


c 

1 ST MAX+KIN 

300010 


XMAX=-1.E10 

OTTO 01 2 


XMTN=i «£1C 

300014 


DO 10 1*1, L 

TUirOiS ' 


TF(2(I> .GT.XMAX) XMAX=Z(I1 

0 COO 22 


IF(Y(I1 .GT.XMAX1 XMAX=Y(I1 

000027 


IF (Z(I) .LT.XMIN) XMIN=7(I> 

C 0 0034 


IF(YCT)«LT* XMIN1 XMIN=Y(I) 

300041 


IF(X(I> .GT.XMAX1 XHAX=X (I) 

0 0 0046 

10 

IF CX(I). LT.XMIN) XMIN*X (1 1 

Ut3 JOHb — 


“ SLC=65./ (XHAX-XHTN) 

00 0 060 


B=-65.*XMAX/(XHAX-XMIN1 *66. 


c 

THEN SCALE 

000064 


XINT* (XMAX-XMIN1/65 . 

000067 


PRINT 99,XMAX,XMIN,XINT 

000101 

99 

FORMAT (IX, *MAX= *,Gi4.6,* MIN* *,G14.6, 



-+* TNT* *, G14.6T 

300101 


IF (XMAX.GT.O..ANO.XMIN.LT.O.) IYZ=INT(B1 

0 00114 


STEEP=0. 

000115 

24 

DO 20 1*1,66 

oooiir 

20 

OUT (I )*STAR 

000123 


IF (XMAX.GT. 0. . AND. XMN.LT. 0.) CUT (IYZI*AZ 

D00IT4 


" IFTSTEER. EC .1 .) GO TO 25 

000136 


OUT (INT (SLO*X (1 1 +B) )*XX 

000144 


OUT (INT (SLO*Y (11 +B1 1 = YY 

000150 


OUT (INT(SLC*Z(11*81 )*ZZ 

0 OH 154 


LL *1 

000155 


PRINT 22,IL,DUT 

FCOIo4 

22 

“FORMAT (IX ,T3» 66 All 

300164 


L ENO*L - 1 

000166 


DO 23 LL=2,LENC 

a o o it 2 


DO 28 1=1,66 

0 001^3 - 

2 8 

OUT II1=8LANK 

000177 


OUT (1I*ST AR 

WOTOO 


0UT(&6T^mn* _ 


TO a 23 0 
300212 
000221 
G0 021JT 
000237 23 

OW2T3 

000 255 

00375 * Z5 

360264 


00 0 273 

0 0 0 30 2 


0 00312 


IF (XMAX.GT.O.. flNO.XHIN.IT. 0.1 OUT (I YZ) =ST AR 

OUT (INT(SLO*X(LL) ♦BIMXX 

OUT (INT (SLO*V (LL)+B))=YY 

OUT (m<SLC*Z (LLWB)) = ZZ 

PRINT 22,LL,OUT 

"3TEFW*T. - 

GO TO 24 

OUT (INT( SL0*X <L) fB) )=XX 
OUT (INT (SLO*V (L1+B) 1=YY 
OUT (I NT (SLO*Z (Li fB) > = ZZ 
PRINT 22, L, OUT 

RTTURU - 

END 


VI 1-2 



0 C 0 006 
' QTTOW6 “ 
OCOOIO 

— 

SUBPOLTINE SURV(CX,V,VS) 
DIMENSION VSC61 
IF (CX .EQ.1.1GO TO 1 

VS (1) =VS (1) +1 . 

0 0 OUii 


T7ST2T~= VS ( 2) +V 

000013 


VS ( 4) = VS (4) +ftBS (V) 

Terrors' 


VS (6)-VS{6)+V*V 

303020 


GO TO 2 

“TXTT0T0 

1 

CONTINUE 

000020 


IFCVSCll.EQ.O.) GOTO 3 

— (Town — 


V5(2)=VS(2V/VS(1) 

300023 


VS (4) =VS (4) /VSC1) 

OUT) 0 25 


0UH=VS<6)/VS<1)-VS<2)*VS 

300030 


IFCOUF.LT .0.1 GCT04 

3 0 3072 


VS(3)=SQRTC0UM) 

300036 


DUN-VS (61/VSC1) -VS (4) *VS 

- ' innnnr? 


X F t UET J • CT *0*1 GCT05 

300044 


VS(5)=SQRT (OUW) 

'0*0 005 0 

... 

PRINT 20, VS 

3000^7 

20 

FORMAT <6E12 *41 

30 30F7' 

" 

RETURN 

003060 

3 

PRINT 10 

' " J 0 0 06TT 

ITT 

FORMAT (1X,*MISSED COUNT* 

3 C 0364 


RETURN 

3 0 306TT 

4 

PRINT 11 

OflOO^l 

11 

FORMAT (IX,* SORT PBLM 1*3 

~~0 0 037T - 

“ 

RETURN 

000072 
"mnnTTK 

5 

“T7 

PPINT 12 

-FTTRMBTClXr^SORT PBLM 2*1 


3 0 00 7 6 2 RETURN 

3 0 OTPT EMU 



SU8R0UTINE CHM£W<S0» 

003003 DIMENSION SB(«,8> 

QOOOOT" **0 . 

000004 AUt*0 * 

TTT3TO5 RUN-fl i 

C C 0006 ALL-0. 
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1 .7*0 Dt + UZ I « ODO0E+O 1 3. *2*4 3E + G z 1*000 OF + 00 1~70 C C O E + Ol — 57090 OE+0O 

FBANP, F«PHA,$0AMP,SOPHA,Gat14 3 f N43, NSLICE 

2 * 0 0 00E+ 00 2.0000E-Q1 2. 30006-01 4.7OQ0E-O1 9. OOOOE-Q2 64 26 

FI 0,FRQOOP,FHDI*T,ASOFT, 11,12, 111,122 

1 ■ OO00E+ 01 0. 0. 2. 64966-01 1111 

SVRHZ,5NHX,$N*N,_ TL 

07 2 . 815OE+04 “9*950 OE+0 3 ~ *8. 0000 E+0 0 

CCl OSC,I DEBIT, GYNR 

0* 1 00 3.0000E+00 

WH,K,TAU1,TAU2^^^ 

1.659344E+02 4.666449E+02 4.499999E-01 8 .442020E-0 3 * 

NNB, K9,T AUB1 , T AU9? '~ - 

9 *42 8 1 6 -01 6.14 8 9 E + 0 2 2 * 0 3536+ 0 2 1 14 9 43 E+0 1 

IBIT,TSYNC 

16 1*13646-02 

PE, 5T0EY,PS5UM t HSSIG,K CUNT 


8.3t58E-t)2 6.33VEE-03 

LOCKIN 2 DROPOUTS 

tPftM 

V.3?5*E-01 2.3<*06E-01 

1AVE 7.^ 00 HE to iso 0. 

2000 






L n nrr 

5.7000E+01 

2.00t0E+tl0 

J.OOOOE+OO 

I.OOOOE*OC 

0. 
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4 * 0 OOOE + 0O 
2*0 0006+00 
1 *00006+00 
0. 

H 

6* 00 0 0E + CO 
4* 000 OE + OO 
0* 

1. 00006+00 
0* 

9. 0000 E+ 00 

0. 

3* 0 Q 00 E+ 00 
0. 

0* 

1. 00006+01 
2.00006+00 
1*00006+00 
0. 

0. 

4* 000 OE + OO 
1 . 00 0 OE + OO 
0. 

6. 

0. 

2 . OOQOE +00 
3.0000E+00 
1.0Q00E+00 
1*000 OE+OO 
0 * 

5*00006+00 

2*00006+00 

0. 

0* 

0* 

2.0000E+00 
1*0 00 OE + OO 

1.00 0 OE+flO 
9* 

1 . 0 0Q0E+Q0 

2 • 00 00E + 00 
3.9O00E+0B 
3* OOOOE+OO 
0* 

0. 

1 • OU UUL* UU 

u • 


~C. 

2.0000F+00 

~0* 

1*00006+00 

0 • 

—07 

9. 

3*9000E+00 

0* 

0* 

0* 
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1 *00006+00 

0* 

0. 

0. 

0* 

0* 

0* 

0. 

0* 

0. 

0* 

0* 

9* 

0. 

0. 

0. 
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0* 

0* 

G* 

0* 

0* 

0. 

0. 

0. 

~0. 

0* 

0* 

0* 

0* 

0. 

0. 

0. 
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639 003 934 1229 

REC SLICED 2.45666+01 3.93886+GC 

XXXC,A5Ut1i, ASUP2, Y5UP1, Y<Ufi2 ,t^UM3 , 

7 *600 2E + 04 1.9156E+06 5.673CE+07 *1 

v h istogra m — 

1 ♦ 4996E-04 
6.8731E-04 
1.2122E-03 
2.3668E-03 
4* 69 07E-O 3 
8 ; 222 7 g - 03 
1.84C2E-Q2 
2*30436-02 
2.6243E-02 
1 *50466-02 


1.4996E-J4 
1.07476-03 
2*03696-03 
3*07426-0 3 
5.3485E-03 
“03 
1*81026- 02 
2.3731E-02 
2*62556-02 
1 *34716-02 


1.4996E-C4 
8.6226E-C4 
1.6121E-03 
3* 1 24 IE- 03 
6*4607 E- 03 
9732236>03^ 
2.1094E-C2 
2. 34 94 E- C2 
2* 79556- 02 
1.09T2E-02 


2 .40 84E+0 7 -4.58856+08 


• 313 3E +06 


1.37466-04 
1 • 1 072E-C3 
1.9745E-C3 
2.9367E-C3 
6*66076-03 
9.62246-63 
1.9232E-C2 
2. 2656E-C2 
2 • 64 30 E- 0 2 
6 • 8351 E- 03 


1.2497E-04 
9.3724E-04 
2*3 1 19E-03 
2 * 3743E— 03 
0.O228E-O3 
1*0 3 35 6-02 

1 *9770 E-Q 2 
2.1944E-02 

2 • 2 219E- 02 
8 * 38526—03 


2 .7492E- 04 
7. 62296-04 
2« 22446-03 
2.72436-03 
7 • 1 655E-03 
1» 073 5E- 0 2 
2.3469E-02 
2.4668E-02 
2 • 0569E-02 
6.9231E-03 


3* 8739E-04 
6.7481E-04 
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4* 2863E-0 3 
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6.4982E-04 
2* 1119E- 03 
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6.3732E-03 
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1.7520E-02 
“ “1*7955E- 03 


7.4979E-04 
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3.4116E-03 
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1*53716-03 5.9904E-O3 
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o. too 3 .ooooe+oo 

MH,K,TIU1,TAU2 _ ~ 
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c » 
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G. 
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0# 
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0. 

0. 
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0. 

0. 

0. 
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0. 

0. 

0. 

0. 

0. 

0. 
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0. 
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0. 
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1. 00006+09 
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0. 
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0. 
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0. 

0. 

0. 

0. 

0. 

0. 

0. 

0, 

0. 

1.90006*00 

0. 

0. 

0* 

0. 

0. 

0. 
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0. 

1.00006*00 
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0. 

9. 

0. 

0. 

1 . 00006+00 

0. 

0. 

0. 
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0. 

1.000 OE +90 

1* 0 OOOE + 0O 
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0. 
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371 87* 1326 *93 

REC St ICE- 2.6266E»01 3.5789E*Q0 

* 7 y C , A SUH 1 , »SU t*2 , V SU t* 1 , ¥SUK? , ¥SU M3 , 
7.600*E*0* 1.915QE*06 5.53S2E*07 - 

1.7*1*E*Q6 

*.1Q5*E*07 -9.885QE*08 





— v — rniaw pw 
1.2fcfl6£-05 
2.2*7*E-0* 
1.373*E-0* 
3.4960E-04 

6.naoe-o* 

1.2486E-fl5 
*2.7468E-04 
1 .87286-04 
4.24516-04 
5 . 3688E-04 

1.2406E-C5 
2.62206-04 
2. 74606- 04 
3.62006-04 
9.61396-04 

■1. 24665-05 
1.99776-04 
2.74606-04 
£.61058-04 
1.12376-03 

1. 24866-05 
3.49606-04 
1.87286-04 
3.496OE-04 
1.01136-03 

1.24066-05 

3.49606-04 

3.24636-04 

5.61856-04 

1.06636-03 

1.2486E-05 
2.6220E-0* 
3.9954E-04 
*• 8694E-0* 
9.3642E-0* 

■"1.2486E-05 
1.9977E-0* 
1.4983E-0* 
5.4937E-0* 
1. 086 3E- 03 

1.24066-04 
1.9977E-04 
3. 7457E-04 
6.11006-04 
1.17 37E-03 

1.2406E-04 

1.0720E-O4 

3.12146-04 

0.30546*04 

1.57326-03 

1.29856-03 

2.64*06-03 
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2.70316-02 

3.08*0E-03 2.73**£-C3 

7.7161E-03 1 • 00 01E-C2 

1.7255E-02 1. B217E-02 
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1.1237E-02 
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2.7955E-02 

“2.20496-03 
3.94556-0 3 
1.12626-02 
2.12636-02 
2.97666-02 

2. 1226E-0 3 
4.0569E-O3 
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2.3710E-02 
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1 * 03546 - 0 3 
5.24406- C3 
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2.4697E-02 
3.25506-02 

1.79796*03 

5.4100E-O3 

1.43036*02 

2.3296E-02 
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1.55956*02 
2.50456-02 
3.33376-01 
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62.0 ,9IF= 1000 ,BL2= 176 ,NSIICE= 26 , ASOFT® .3000 

= 2.ai S 0 F» 04 t SMWM a 6.8 S 00E + 03 , F LO= l.Q O OOE +O l , FR( 

,SOAMP= .230 ,SDPHA= .470 ,11,12,111,122 129 130 


3.6842E-02 ,STDEV® 4.3216E-03 ,¥ MEAN® 4.6193E-01V STO OEV= 2.0616E-01 ,B1TS= 1900 

1 r ■» /I ^ nn. , n«rnr T ur 1 ,1 i f nnnr * (i , T 


LOCKING i 9 DROPOUTS 3 1 ,AVE= 2.8000E401 , STO DEV= 0. 


3.100QE>01 4 • Q 0 00 £♦ 0 0 3. 0000E+00 2.0000E400 2.0000 E40Q l.QOQOE+OD 


2.0000E+00 


1.0OQOE+OO 


!• QQ00E + 00 0« 

















E NO 08= 9.4 »9£LF= 62.0 ,BIF= 1000 ,BL2= 176 ,NSLICE = 26 ,ASOFT= .3000 

, SW 4 K= — 2« 0 1 S 0E * B 4 ,S WMN = — 6. B 5 00E + 0 3 , FLO= — l.Q000E » 01 , FRQDOP 

»OYNR= 9.0 ,SDAt1P= .230 t SOPH A= .470 ♦11,12,111,122 129 130 

n a £. T lilt — l C ft nrt r rtj a i. ^ ft r n -r uuO - n i/n . r 4 n 


PE= 2.4737E-J2 f STDEV= 3.5633E-03 f V MEAK= 4* 56 33E-I3L V STQ DEV= 2.0596E-G1 ,BITS 


LOCKING 1 t DROPOUTS- i « AVE= 4.0000E+00 , STD 0EV= 0* 

-E*** 

8.0000E*00 3.000QE»OQ 1.0050E*00 1.0000E+00 1.003 CE+0Q 0* 3* D 00 OE+QG 1# QQ0QE+0Q 1« 000QE+0O 

— i *o a ofl € >oq — i .aaao € *oo — — i^oooo e» oo — < u — — < u 

0. 1 « OGOQE+QQ 0. 0. 0. l*OOQOE*O0 0. 0* 0* 
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a. — — 

1* OQG0E+OO Oo 
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ENODB* 9.4 »OELF* 62.0 ,BIF= 1000 ,BL2= 176 ,NSLICE= 26 , ASOFT* .3000 


TL- -6.50 ,DVNR= 


|^^||| 


SDAXP* .230 , SOPHA* .470 ,11,12,111,122 129 130 


PE= 2.8421E-02 ,STO£V* 3.S123E-03 ,V HEAR* 4.6307E-01V STD DEV* 2.0835E-Q1 ,BITS 


LOCKIN* 0 , DROPOUTS* 1 ,AVE= S.OOOOEtOO , STD OEV* 0. 



1.40Q0E+01 3.00 00£*00 3.0000E*00 0. 


0. 1.0000E*00 l.OOQOE+OO 1.000QE+00 

1. 00 ODE* 00 1.0000E*00 0. 0. 


1.0000E+00 0. 


2.00006*00 
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ENQOM 9.4 ,0ELF= 62.0 ,BIF= 1000 ,BL2= 176 ,NSLICE= 25 , ASOFT= .3000 

SVR H Z- 0 . ,S H MK« - 2 .6 1 50 E + 0 4 , SW W » — 6 . 85 00E + 0S , FLO» 1. ft OBOE Ml * FRQOQP= — 8 . ,F ft O P RT= l ,0ft 00E »0 1 

TL= -6.50 , OVNRs 0.0 , SQAMP= 0.000 ,SOPHA= 0.000 ,11,12,111,122 129 130 1 1 

-WN - 1 65,9 , «■ 468 . 6 , T A01« — 4.50 00E - 01 ,TAU2» — 8.44 2 0 E -83 ,HM9~ .9 »K8= 6 1 4 .9 , TAU 8 1= 2.03 5 3E + 02 , T A U ft 2s — 1 .4 9 4 SE MO - 
DE* 1.0526E-02 ,ST0€V* 2.3413E-03 ,\f MEAN* 4.7068E-QLV STO DE¥= 1.9479E-01 »BITS= 1900 

-P Ef r 6T - 0 6ft» — 1 .8 17 4 E - 0 3 , I H £ AM = I STO 06V» — 5-. 8« 07E -Mft-r»Z^S- - — ft. 

LOCKIN= 1 ,9ROPOUTS=0 

ERR M — — — 

1.0000E+00 0. 0. 0. 0. 0. 0. 0. 0. ft. 

ft^- -ft, ft. - ft, l .ftftftB EM ft ft, — 8. 8. 0, 8. 

8. 0. 0. 1.0008E+00 1.0000E+80 0. 0. 0. 0. 1.OOO0E+OO 

ft, — — B-, ft, ft. 8. ft, ft. ft- ft. ft. 

9. a. 1.000QE+Q8 0. D. 0. 1.0000E+00 0. 0. 1.O0OOE+0O 

9, 1 . 0 8 0 8 E +88-2 . B 00 8 E + 8 8- ft, ft. ft. ft. 8. ft. »• 

0. 0. 0. 0. 0. 0. 0. 1.000QE+00 0. 0. 

9, 4, --ft, ft. ft, ft. 8 , ft. ft. ; 8. 

ft. 0. 0. 0. 0. 0. 0. 0. 8. ft* 

ft, ft, ft, ft. —ft. ft. ft, ft. ft, 8, 00 08 € »0ft 


— ftr-1 Pt 8 / A >TI H E 5 451 

lftZ 26 9 4 0 1 1 

9ft 3« tft 6 3 ft- ft— 

42 30 26 8 1 0 0 

1-0 1ft 9 3 1- — 0 8 

2 0 2 2 0 0 8 

4 4 ft ft 4 - —ft ft 

0 ft 0.0 0 0 0 

ft. 0 ft 4 4 ft - — O — 

+1,-1 P(3/A)TIMES 482 

1 1 0 3 16 6ft 125 — 

0 0 0 9 2 11 20 

ft-- —ft 8 0 2 4 - 4 — 

0 0 0 0 0 0 3 

4 B- ft 0 ft 1 0 — 

0 0 0 0 1 0 0 

8 0 0 9 1 8 0 — 

0 0 0 0 8 0 0 

11 , 1 2, I ll,I2g 1594 — 1314 1 1 

REC SL ICE- 2.5830E+01 ,STD DEV= 3.5020E+00 
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0- 

0 

— 0- 

0 
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46ft- 
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12 

3 

0 

0 
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ft ft— ft- 9 ft ft 0 ft 
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1 q 1 ft o 0 0 0 
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ftft 7 4 2 -ft — -ft —ft 0 

69 29 9 3 1 0 0 0 

149 46 53 15 1- ft ft 0 

♦1,+1 P (B/ A) TINES 485 

0 ft ft 8 8 2 0 0 

00000000 

4 ft 0 ..... 0 o — — 1 0 ft 

0 0 0 0 0 0 0 2 

0 0 0 ft ft 3 7 10 

0 0 0 0 5 6 31 60 

0 e 4 e 3 39 52 ftft 

1 1 1 1 4 25 41 111 

■ - DYNR = 0 

2/7/74N 
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APPENDIX IX 


PARAMETER VARIATION PRINTOUTS 


TITLE PAGE 

DELF = 26.4, Fading IX-2 

DELF = 44, " 3 

DELF = 88 , 4 

DELF = 132, " 5 

DELF = 26.4, No Fade 6 

DELF = 44, 7 

DELF = 88 , 8 

DELF = 132, " 9 

B1F = 500, Fading 10 

B1F = 750, " 11 

B1F = 1500, " 12 

B1F = 2000, " 13 

B1F = 500, No Fading 14 

B1F = 750, " 15 

B1 F = 1500, " 16 

B1F = 2000, " 17 

BL2 = 88 , Fading 18 

BL2 = 100, " 19 

BL2 = 164, " 20 

BL2 = 188, " 21 

BL2 = 88 , No Fade 22 

BL2 = 100, " 23 

BL2 = 164, " 24 

BL 2 = 188, % 25 

TAU4 = 1.8 x 10,, Fading 26 

TAU4 = 1.1 x 10"f, " 27 

TAU4 = 7.9 x 10",, " 28 

TAU4 = 1.8 x 10",, No Fading 29 

TAU4 =1.1 x 10 4 , 11 30 

TAU4 - 7.9 x 10 , 11 31 

WORST CASE 1, Fading 32 

WORST CASE 2, 11 33 

WORST CASE 1, No Fading 34 

WORST CASE 2 , 35 
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TL= -6.50 ,0YMR= 3.0 ,SDA*P = .230 , 
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ft 
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■ — • 

0 o 
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45 

53 

27 
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1 
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ez 
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7 
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3 
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21 

4 
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0 

* 

z 
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0 
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4 
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- - 0 

- 1 

D 0 
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5.0QOOE + ftft 
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4 oOoooe+eft 
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APPENDIX X 


ALTERNATE CONTINUOUS PHASE MODULATORS 

Three alternate continuous phase modulators, a three step upconverter, 
an indirect FM modulator and a single step upconverter were considered during 
the course of this study. The latter was finally selected because of its 
simplicity and minimum parts count. This appendix discusses these alternate 
designs. 

Three Step Upconverter 

Figure 1 shows the block diagram for the modulator and upconverter. The 
400 MHz FSK signal is synthesized by generating a 1 MHz signal in a voltage 
controlled crystal oscillator and up-converting this signal by mixing it with 
higher frequencies derived from a 23 MHz crystal oscillator. The mixing is done 
in three steps to make filtering practical at each frequency. FSK is accom- 
plished by pulling the 1 MHz VCXO with a +3.2 volt DC signal amplified from the 
TTL data. 

The oscillators are specified to have an end-of-life long term drift of 
+10 ppm. The 1 MHz oscillator is voltage controlled to provide the 64 Hz 
frequency shift. A VCO was chosen over frequency switching to provide phase 
continuity between frequencies. 

A typical modulation range for a crystal oscillator is 80 ppm (Frequency 
Electronics, Inc.). For the needed 64 Hz modulation, this requires a 1 MHz 
center frequency and a DC control voltage of +3.2 volts. The frequency will 
typically switch in 20 microseconds. 

The power losses and gains are indicated in Figure X-l. The criteria for 
the design are: 

o An input of 6 dBm into the local oscillator port of the mixer. 

o 0 dBm or less into the signal port to maintain at least a 6 dB ratio 
to minimize spurs. 

o Limit amplifier gain to approximately 10 dB due to environmental 
constraints. 

The mixers have a loss of approximately 6 dB at a ratio of 6 dB between 
signal and L.O. (per Relcom Co.). 


X-l 



Filters are provided in the L.O. lines to eliminate the harmonics of the 
square waves generated by multiplying and dividing. This application is not 
critical . 

Filters are also provided after each mixing operation to eliminate harmonic 
spurs. The filter bandwidths were chosen to have a -60 dB bandwidth at less 
than the separation of the closest low order spurs. 

The spurs of the first mixer were the closest and a crystal filter is 
required. The crystal filter is specified by Crystal Network Products. The 
2nd and 3rd mixer filters are designed according to Texscan tubular filter 
specifications and are summarized in Figure 1. The main objection to the 
Figure 1 design* is its overall hardware complexity, high parts count, and con- 
sequently, its lower reliability, higher power consumption, and added weight. 
Based on information obtained from a crystal oscillator manufacturer alternate 
designs are possible. 

The following are the constraints which dictated the design parameters: 

1. The modulation index is .7 radians +5% over the temperature range 
of -40°F to +160°F. 

2. Over the temperature range, an ovenized oscillator can be pulled a 
minimum of +1 ppm and a maximum of +100 ppm. The minimum value is 
indicated by a crystal oscillator vendor as being the value one could 
pull an oscillator to an accuracy of <5%. The maximum is constrained 
by the phase noise of an oscillator which should be >76 dBc at 100 Hz 
from the carrier in a 1 Hz bandwidth. To deviate +31 Hz, this implies 
a VCXO may be centered between 300 KHz and 30 MHz. However, due to 
the limitations of practical crystal cuts, the best oscillators are 
constructed in the 5 to 20 MHz range. 

3. The minimum two-sided bandpass filter 3 dB bandwidth is to be 3%. 

The filter attenuation of the closest spur should be at least 4-0 dB 
down at the modulator output. This is because the non-linear power 
amplifiers will enhance the spurs approximately 15 dB and also create 
intermodulation products. Because of the extreme temperature range 
and associated filter drift problems, 60 dB will be used as the design 


*M. D. Saferstain, "SAEP Modulator", TRW IOC 7322.4-15, November 19, 1973. 
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constraint. The filter will have a maximum of 5 poles to minimize 
filter size and number of elements and maximize reliability. 

Given these constraints, two alternate modulator configurations are real- 
izable as shown in Figures X-2 and X-3. 

Indirect FM Modulator 

The indirect FM modulator design based on the narrowband phase modulator 
concept conceived by Armstrong, integrates the data, linearly phase modulates 
the resultant waveform with a modulation index of 0.0775, and multiplies the 
modulated carrier 8 times to achieve an output frequency of 400 MHz at a modu- 
lation index of 0.7. Linear phase modulation is achieved by summing the 
carrier in quadrature with the double sideband suppressed carrier modulated 
signal; i.e., the output Y(t) as a function of the input signal X(t) is, 

Y(t) = cos (w c t) + X(t)sin(w c t) 

= cos(u) c t + tan _1 X(t)) - cos(a) c t + X(t)) 
for jX(t)| «1. In this case, X ( t ) = 3/<j>(t)dt where 3 is the modulation index 



FIGURE X-l 

THREE STEP UPCONVERTER MODULATOR 
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LOCAL : 

OSCILLATOR 
50 MHz 

FIGURE X-2 

INDIRECT FM MODULATOR 

and <|>(t) is the digital data. The one fundamental drawback to this concept is 
the practical implementation of the integrator and mixer. For long strings of 
marks or spaces, or where the average number of marks and spaces are not equal, 
the integrator and mixer must have sufficient dynamic range to operate over the 
entire period of data transmission. Furthermore, the phase modulation^apprgxi- 
mation holds only for the range where | X ( t ) [<< 1 . Since tan~^a = a - . . 

if we modulate one bit at a modulation index of .0775, a string of data where 
the marks outnumber the spaces by 10 bits would increase the distortion from 
.2% to 20%**, clearly an unacceptable situation. Lower modulation indices may 
be possible by decreasing the 50 MHz L0 and increasing the multiplication fac- 
tor. However, in view of the simplicity of the single conversion direct FM 
design discussed below, this design is untenable, 

**"The desired signal at the output of an FM modulator is: 

Y(t) = cos(w ct + e/4>(t)dt). 

The actual signal out of the narrowband FM modulator is: 

Y (t) = cos [w t + tan -1 (e/4>{t)dt)] 

L 35 

= cos ^ e/4> ( t)dt - ^-3- {/<f)(t)dt)3 + ( )** -+***] 

where the higher power terms in the expansion are the distortion terms. If for 
one bit, e/$(t)dt = .0775 radians, then the distortion is in the order of 
(.0775) /3 = .2%. Where the data is such that / <f>(t)dt = 10 /<t>(t)dt, the dis- 
tortion increases to ( . 775 ) 2 /3 * 20%. 
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Single Conversion FM Modulator 

The single conversion direct FM. modulator shown in Figure 3 takes the 0 
dBm output of a FCXO at 20 MHz and upcon verts in one step to the output 
frequency of 400 MHz. The nearest mixer sideband in the upconversion is 40 MHz 
away using a doubly balanced mixer which also suppresses the 380 MHz L0 a 
minimum of 25 dB. This implies a 3% bandwidth, 5 pole Chebyshev filter may 
be used at 400 MHz which would result in the 380 MHz L0 being attenuated an 
additional 40 dB to a total level of 65 dB and the 40 MHz sideband attenuated 
to 68 dB. The total loss in the filter is approximately 5 dB, thus an 11 dB 
gain amplifier is needed to bring the output to 0 dBm. This modulation techni- 
que is simple and straightforward, requiring a minimum of parts, meets all the 
design constraints and, based on the new oscillator data, is the TRW recommended 
design. 
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FIGURE X-3 

SINGLE CONVERSION DIRECT FM MODULATOR 




